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Abstract 

We present a review of scattering models that can be used to describe the low-energy 
behavior of identical bosonic atoms. In the simplest models, the only degrees of free- 
dom are atoms in the same spin state. More elaborate models have other degrees of 
freedom, such as atoms in other spin states or diatomic molecules. The parameters 
of the scattering models are specified by giving the S-wave phase shifts for scatter- 
ing of atoms in the spin state of primary interest. The models are formulated as 
local quantum field theories and the renormalization of their coupling constants is 
determined. Some of the parameters can be constrained by renormalizability or by 
the absence of negative-norm states. The Green's functions that describe the evo- 
lution of 2-atom states are determined analytically. They are used to determine the 
T-matrix elements for atom-atom scattering and the binding energies of diatomic 
molecules. The scattering models all exhibit universal behavior as the scattering 
length in a specific spin state becomes large. 



1 Introduction 



The development of methods for trapping atoms and cooling them to ultralow 
temperatures has launched the new field of cold atom physics. Among the 
dramatic achievements in this field have been the creation of Bose-Einstein 
condensates of bosonic atoms and superfluids of fermionic atoms. At the low 
temperatures at which these phenomena occur, the de Broglie wavelengths of 
the atoms are much larger than their sizes or the ranges of their interactions. 
This makes the behavior of the atoms insensitive to many of the details of 
atomic structure and interatomic interactions. The atoms can be treated as 
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point particles and their interactions can be described very accurately by 
simple models. 

One class of models that is particularly convenient for describing ultracold 
atoms is scattering models. Scattering models are formulated in terms of pa- 
rameters that describe the low-energy scattering of atoms. The simplest exam- 
ple is the Zero-Range Model whose only interaction parameter in the 2-body 
sector is the S-wave scattering length. The scattering length is the most impor- 
tant interaction parameter governing the behavior of cold atoms. The primary 
advantage of scattering models for describing cold atoms is that their param- 
eters are directly related to physical observables at the relevant energy scale. 

Another advantage of scattering models is that there are many useful mod- 
els that can be solved analytically in the 2-body sector. There are very few 
analytic results on the 3-body and higher n-body problem, but having an ana- 
lytic solution to the 2-body problem is a great simplification in the numerical 
solution of higher few-body problems. The analytic solution of the 2-body 
problem is also useful in the many-body problem, because interaction effects 
in low-density systems of atoms are dominated by 2-body interactions. 

A particularly convenient class of scattering models is ones with local inter- 
actions. In these models, atoms interact with each other only when they are 
at the same point in space. Such a simplified representation of the interac- 
tions for cold atoms is useful because the large de Broglie wavelengths of the 
atoms makes them insensitive to the range of interatomic interactions. As a 
consequence, the low-energy behavior of the atoms can be described to a good 
approximation by zero-range interactions. 

A scattering model with local interactions can be described by a local quantum 
field theory. A formulation in terms of a quantum field theory is useful because 
the same framework describes the 2-body, higher few-body, and many-body 
problems. Quantum field theories with local interaction terms are particularly 
convenient, because theoretical methods for dealing with such theories are 
well developed. Much of the stimulus for this development has come from the 
success of local relativistic quantum field theories in describing elementary 
particles. One of the complications of local quantum field theories is that they 
are inherently singular at short distances, but this problem can be handled 
using the machinery of renormalization. 

In this review, we present a unified treatment of scattering models that can 
be used to describe the low-energy behavior of atoms or other nonrelativistic 
particles with short-range interactions. We restrict our attention for the most 
part to bosons with the same mass and with S-wave interactions only. The 
extensions to fermions, to particles with different masses, and to interactions 
in higher angular momentum channels are straightforward using the formalism 



2 



of quantum field theory. In the simplest scattering models, the only degrees 
of freedom are atoms that are all in the same spin state. We also consider 
scattering models whose degrees of freedom include atoms in other spin states 
or diatomic molecules. 

We define the physical parameters of each scattering model by specifying the S- 
wave phase shift for the spin state of primary interest. We formulate the model 
as a local quantum field theory with an ultraviolet cutoff. The renormalizations 
that relate the coupling constants of the quantum field theory to the physical 
parameters of the scattering model are determined. The Green's function that 
describes the evolution of a 2-atom system is derived analytically. It is used 
to determine the T-matrix elements for atom-atom scattering and the binding 
energies of diatomic molecules. 

We begin the review in Sec. |2]by describing a fundamental theory that provides 
an extremely accurate description of the low-energy behavior of alkali atoms. 
We then give a general discussion of scattering models in Sec. El In Sec. HI we 
discuss the Zero- Range Model, whose only interaction parameter in the 2-atom 
sector is the scattering length. In Sec. El we discuss the Effective Range Model, 
which has a second interaction parameter. Other scattering models whose only 
degrees of freedom are an atom in a single spin state are described briefly in 
Sec. El In Sec. El we discuss the Two- Channel Model, which describes atoms 
that can be in either of two spin states. In Sec. El we describe the Resonance 
Model, in which a diatomic molecule enters as an additional degree of freedom. 



2 The Fundamental Theory 

In this section, we describe the fundamental theory that provides an extremely 
accurate description of atoms with energies small compared to the splitting 
between the ground state of the atom and its first electronic excitation. To be 
specific, we focus on alkali atoms. These atoms have been used extensively in 
cold atom experiments, because they have properties that make it particularly 
easy to cool them to ultralow temperatures using current technology, such as 
laser cooling and evaporative cooling. 

2.1 Hamiltonian 

The alkali atoms H, Li, Na, K, Rb, Cs, and Fr can be labeled by an integer 
n = 1, 2, 3, 4, 5, 6, and 7. The electronic structure of the n th alkali atom 
can be roughly approximated by a single valence electron in the nS state of 
a Coulomb field created by the closed shells of the inner electrons and the 
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nucleus, which together have total electric charge +1. Thus the energy scale 
of the first electronic excitation of the Rb atom can be estimated from the 
hydrogen spectrum: Electronic ~ E Rydhcrg /n 2 , where -E Ryd berg = 13.6 eV. 

An alkali atom in its electronic ground state has multiple spin states. There 
are two contributions to its spin: the electronic spin S with quantum number 
s = | and the nuclear spin J with quantum number i. The 2(2i + 1) spin 
states can be labeled \m s ,mi), where m s and mi specify the eigenvalues of S z 
and I z . The Hamiltonian for a single atom includes a hyperfine term that can 
be expressed in the form 



2E M 

^hyperfine ^^^f^ 1 ^- W 

This term splits the ground state of the atom into two hyperfine multiplets 
with energies differing by E^. The eigenstates can be labeled by the eigenval- 
ues of the hyperfine spin F = I + S. The associated quantum numbers / and 
rrif specify the eigenvalues of F 2 and F z . The eigenvalues of -ffhyperfine are 



E f,m f ~ M - [ ) 

The hyperfine multiplet consists of 2i + 2 states with / = i + | and 2i states 
with / — i — |. 

In the presence of a magnetic field B = Bz, the Hamiltonian for a single atom 
has a magnetic term. The magnetic moment n of the atom is dominated by 
the term proportional to the spin of the electron: /j, — /i S/(\h). The magnetic 
term in the Hamiltonian can be expressed in the form 



-^magnetic — % ^ ' ^ 

If B ^ 0, this term splits the two hyperfine multiplets into 2(2z + 1) hyperfine 
states. In a weak magnetic field satisfying fiB <^ Em, each hyperfine multiplet 
is split into 2f + 1 equally-spaced Zeeman levels \f,rrif). In a strong magnetic 
field satisfying fiB ^> E^, the states are split into a set of 2i + 1 states with 
m s = — i whose energies increase linearly with B and a set of 2i + 1 states 
with m s = +2 whose energies decrease linearly with B. Each of those states 
is the continuation in B of a specific hyperfine state \f,m,f) at small B. It is 
convenient to label the states by the hyperfine quantum number / and mj 
for general B, in spite of the fact that those states are not eigenstates of F 2 
if B 7^ 0. We denote the eigenstates of -f^hypcrfme + -^magnetic by \f,rrif;B) and 
their eigenvalues by E^ mf {B). The two eigenstates with the maximal value of 
|m/| are independent of B: 
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Fig. 1. The hyperfine energy levels as a function of the magnetic field for an alkali 
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atom with i = |, such as 87 Rb 



/ = i + \,m f = ±(i + l);B) = m s = ±\,mi = ±ij . 
Their eigenvalues are exactly linear in B: 



E f,m f {B) 



2i + l 



E hi ^ \iB. 



(5) 



If B ^ 0, each of the other eigenstates \f,rrif]B) is a linear superposition of 
the two states |/ = i — |, rrif) and |/ = % + |, rrif). 

As an illustration, we take 87 Rb atoms, whose nuclear spin quantum number 
is % — |. Since Rb is the 5 th alkali atom, the energy scale of its first electronic 
excitation is approximately -ER y dber g /5 2 ~ 0.55 eV. At B = 0, the electronic 
ground state is split into hyperfine multiplets with f — 1 and / = 2, with the 
f — 2 multiplet higher in energy by = 2.83 x 10~ 5 eV. Thus the hyperfine 
splitting is more than 4 orders of magnitude smaller than the energy scale for 
electronic excitations. The magnetic moment \i of a Rb atom is approximately 
that of an electron: \x ~ 2/xg, where \x& is the Bohr magneton. The magnetic 
energy scale fiB is comparable to when B is about 2400 Gauss. The 
dependence of the hyperfine energy levels on the magnetic field is illustrated 
in Fig. CD 



If the atoms have sufficiently low energy that none of their electronic excita- 
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tions can be excited, they can be described as point particles with multiple 
spin states that interact nonlocally through a potential. In the case of two 
alkali atoms, their interactions at sufficiently low energies can be described 
by the Born-Oppenheimer potentials V s (r) and V t {r) for atoms whose valence 
electrons are in spin-singlet and spin-triplet states, respectively. The spin- 
dependent potential can be written as an outer-product operator acting on 
the space of 2-atom spin states \m s , mi) ® \m' s , m'j): 



V{r)=V t {r) 



i) + (EiS* ® s^/h 2 



+V s (r) [J(1®1)- (EiS^S^/h 2 
The asymptotic behaviors of the two potentials at large r are 
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V t (r), K(r)— E --|, (7) 

where Eq is the scattering threshold for two atoms in the absence of the 
hyperfine and magnetic interactions. The scattering thresholds are split by the 
hyperfine and magnetic interactions in Eqs. ([T]) and Q. The threshold for two 
atoms in the states \f,mf]B) and \f',my,B) is 2Eq + Ef tTnf (B) + Efi tTn > (B). 

There are also 3-body and higher n-body interactions between atoms. Their 
effects are generally believed to be negligible in most cases, and we will ignore 
them. We will also ignore relativistic effects and retardation effects, which 
change the ultimate asymptotic behavior of the van der Waals potential from 
1/r 6 to 1/r 7 . 

The nucleus of an atom with atomic number Z and atomic mass number A 
contains Z protons and A — Z neutrons. The atom also contains Z electrons to 
neutralize the electric charge. Thus the total number of fermionic constituents 
in the atom is A + Z. The atom is a boson if A + Z is even and a fermion if 
A + Z is odd. In the case of alkali atoms, the atomic number Z is always odd. 
Thus an alkali atom is a boson if A is odd and a fermion if A is even. 

If two atoms are in the same spin state \ f,rrif;B), their wavefunction VK r i> r 2) 
must satisfy a symmetrization condition. If the atom is a boson, the wavefunc- 
tion must be symmetric under interchange of the two coordinates: ip(r2, Ti) = 
+ijj(ri, r 2 ). If the atom is a fermion, the wavefunction must be antisymmetric 
under interchange of the two coordinates. More generally, two atoms can be in 
superpositions of the various hyperfine states, in which case their wavefunction 
has multiple components ip m , al m il ,m s2 m i 2 ^2)- The symmetrization condition 
must be applied to the spin quantum numbers as well as to the coordinate 
dependence of the wavefunction. If the atom is a boson, the components of 
the wavefunction must satisfy 
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' l Pm s 2m i 2,m sl mi 1 m s im,n, m a 2mi2 

We can now specify the fundamental theory. The Hamiltonian is the sum of a 
one-body term for every particle and a 2-body term for every pair of particles. 
The one-body term is the sum of the kinetic energy P 2 /(2m), the hyperfme 
term in Eq. (pQ), and the magnetic term in Eq. ([3]). The 2-body potentials are 
given in Eq. ([6]). The bosonic or fermionic nature of the atoms is implemented 
through constraints on the states. If the atom is a boson, the iV-particle wave- 
function must be symmetric under interchange of any pair of atoms. If the 
atom is a fermion, the iV-particle wavefunction must be antisymmetric under 
such an interchange. 



2.2 Quantum field theory formulation 



The fundamental theory can be formulated as a quantum field theory. This 
formulation involves 2(2* + 1) quantum fields 4>m s mA r ) that can be arranged 
into a column vector ty(r). The Hamiltonian is 



Hf un = [ o?r (— W f • W + -^-$t j . S m - 2fiB ■ 

J \2m 2i + 1 / 

+\ J d 3 n J d 3 r 2 (tf (n) <g> ^(r 2 )) f V(r 12 ) (*(n) ® *(r 2 )) , (9) 

where ri 2 = \ fi — t*2| - We have set ft = 1 to simplify the expression. Dimen- 
sional analysis can be used to reintroduce ft if desired. The quantum fields 
ipm a mA r ) satisfy equal-time commutation relations. If the atom is a boson, 
the commutation relations are 



tpm s mA r )^m' s m' t (r')\ = 0, (10a) 
^m s m i ( r )^L'm'.( r ')} = <*m s m^ mim ^ 3 (r - r'). (10b) 

If the atom is a fermion, the commutation relations are replaced by anti- 
commutation relations. The commutation relations imply that the quantum 
field ip msmi ( r ) annihilates an atom in the spin state \mi,m s ) at the point 
r. They also enforce the constraint that quantum states must be symmetric 
under interchange of two atoms in the same spin state. 

The fundamental theory has several symmetries: 

• phase symmetry. The Hamiltonian H{ un is invariant under phase transforma- 
tions: \&(r) — > e l6 ^>(r). This symmetry is associated with the conservation 
of the number of atoms. 
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• translational symmetry. If the magnetic field B(t) is homogeneous, Hf un is 
invariant under translations in space: r — > r + a. This implies the conser- 
vation of the total momentum P. 

• time translational symmetry. If the magnetic field B(r) is static, the theory 
is invariant under translations in time: t — > t + a. This implies the conser- 
vation of the total energy E, which is the eigenvalue of the Hamiltonian 



• rotational symmetry. The Hamiltonian is invariant under rotations gener- 
ated by the total angular momentum operator J = L + S + I, where L is 
the orbital angular momentum, provided the magnetic field B is also trans- 
formed by the rotation. If B — 0, this symmetry implies the conservation 
of the total angular momentum J. If B = Bz, it implies the conservation 
of the component J z . 

• parity symmetry. The Hamiltonian is invariant under the reflection r — > — r. 

The fundamental theory also has Galilean symmetry. Since the Hamiltonian 
is the generator of translations in time, the column vector of quantum fields 
can be extended to a time-dependent operator ty(r,t): 



The transformation of this time-dependent operator under a Galilean boost 
with velocity vector v is 



This symmetry is not associated with a conservation law, because its genera- 
tors do not commute with the Hamiltonian. 

If we restrict our attention to low energies with respect to a given scattering 
threshold, the fundamental theory has additional approximate symmetries. 
If the atoms have sufficiently low kinetic energy, their total orbital angular 
momentum L can be neglected. If B — 0, the total hyperfine spin F = I + S 
is then conserved. If there is a nonzero magnetic field B = Bz, the component 
F z is conserved. The restriction on the energy may also impose a restriction 
on the possible spin states of the atoms. In this case, the conservation of 
total atom number may be replaced by a more restrictive conservation law. If 
B — 0, the energy restriction could allow only atoms in a specific hyperfine 
multiplet with quantum number /, in which case the total number of atoms 
in those 2f + 1 hyperfine states will be conserved. If B is large, the energy 
restriction may require all the atoms to be in the same spin state, in which 
case the number of atoms in that spin state will be conserved. 




(11) 
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2.3 Scattering 



If two atoms in the hyperfine states \ fi,rrifi] B) and \ f2i m f2] B) collide, they 
can either scatter elastically or they can scatter into a different pair of hyper- 
fine states \ f[, m'fi, B) and \ f 2 , m'f 2 ] B). The scattering rate for each scattering 
channel is determined by the T-matrix element, which is a function of the rel- 
ative wavenumbers k and k! for the initial and final states. If the collision 
energy relative to the scattering threshold is sufficiently small, elastic scatter- 
ing is dominated by the 5-wave orbital angular momentum channel and the 
T-matrix element Tf imfl j 2mf2 (k) for elastic scattering reduces to a function 
of k — \k\ — \k'\ only The low-energy limit of the T-matrix element is given 
by the scattering lengths for the hyperfine channel. If the two hyperfine states 
are distinct, we denote the scattering length by o/ im/11 / am „. The low-energy 
limit of the T-matrix element is 

1~fimf 1 ,f2mf2\k) * ~~ fl /i m /iJ2 m /2 as /c — > 0. (13) 

If the two atoms are in the same hyperfine state \f,rrif;B), we denote the 
scattering length by a f m , . The low-energy limit of the T-matrix element is 

87T 

Tfm f jm f (k) — ► —a fmj as k -> 0. (14) 

If the atoms are fermions, the scattering length a/ m , vanishes because the 
wavefunction for an S-wave state must be symmetric under interchange of the 
coordinates. The cross section for the scattering of distinct hyperfine states are 
obtained by squaring the T-matrix element Tf imfl j 2mf2 (k) and multiplying it 
by the flux factor m/(2k) and by the phase space factor mk/(27r). If the two 
atoms are in the same hyperfine state, the phase space factor is multiplied 
by a factor of | to compensate for the double- counting of the states of the 
identical particles. Thus the low-energy limit of the cross section for the elastic 
scattering of two atoms in the \ f,rrif] B) hyperfine state is 

creia S tic(£) — ► 87r|a /m/ | 2 as E -> 0. (15) 

In the absence of the hyperfine and magnetic interactions, scattering through 
the spin-singlet and spin-triplet potentials in Eq. ([HD determines the spin- 
singlet and spin-triplet scattering lengths a s and a t . The spin-triplet scattering 
length a t is the scattering length for a pair of atoms that are both in the same 
hyperfine state given in Eq. (j3J). The only effect of the hyperfine and magnetic 
terms on this channel is to change the scattering threshold. If it were not for 
the hyperfine term in Eq. 01) , the scattering lengths for all the other hyperfine 
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Fig. 2. Qualitative behavior of the spin-singlet potential (lower solid line) and the 
spin-triplet potential (upper solid line). The low-energy scattering properties can 
be reproduced by replacing both potentials by the van der Waals potential —Cq/t 6 
with a boundary condition at r = ro such that the spin-singlet and spin-triplet 
scattering lengths have the correct values a s and at- 

channels could be expressed as linear combinations of a s and a t using Clebsch- 
Gordon coefficients. Because of the hyperfine term, it is necessary to solve the 
coupled-channel Schrodinger equation numerically to determine the scattering 
lengths for all the other hyperfine channels. However, if the scattering lengths 
a s and a t are known, there is a simple approximation to the coupled-channel 
problem that reproduces the low-energy scattering observables, including the 
scattering lengths, quite accurately [lj . As illustrated in Fig. El the spin-singlet 
and spin-triplet potentials V s {r) and V t {r) in Eq. ([6]) can be replaced by their 
asymptotic forms in Eq. ([7]) together with boundary conditions at a small 
separation r . The boundary conditions are tuned so that the spin-singlet and 
spin-triplet scattering lengths have the desired values a s and a t for B = and 
£m — 0. Those same boundary conditions are then used for nonzero values 
of B and E^. Thus the low-energy scattering problem for all the hyperfine 
channels is essentially determined by just 5 parameters: a s , a t , Cq, Em, and 
fiB. 

We now focus on a specific hyperfine channel. To be definite, we take the alkali 
atom to be a boson. We also take the two atoms to be in the same hyperfine 
state \f,rrif;B), so they are identical bosons. We choose the origin of energy to 
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be the scattering threshold for that channel: 2Eq + 2Ef )m .(B). The T- matrix 
element for S-wave scattering of two particles with collision energy E = k 2 /m 
can be expressed in the form 



where So(k) is the S-wave phase shift. If there are no lower hyperfine channels 
into which the atoms can scatter, the unitarity of the S-matrix implies that 
the phase shift is real valued. The low-energy behavior of the phase shift is 
given by the effective range expansion: 

k cot 5 (k) = -1/a + \r s k 2 + . . . , (17) 

which defines the scattering length a = a,f m , and the effective range r s . Tak- 
ing the limit k — > in Eq. (fl6l) and using Eq. (fTTl) . we recover the limiting 
expression for the T-matrix element in Eq. (TTJ 



If the fundamental theory had only short-range potentials, kcot 8o{k) would 
be an analytic function of k 2 at k = and it would therefore have a power 
series expansion in k 2 . However, the asymptotic behavior of the potentials 
V s (r) and V t (r), as given in Eq. ([7j), is 1/r 6 . As a consequence, beginning at 
order fc 4 , there are logarithms of k in the coefficients of the powers of k. Thus 
k cot So(k) has a well-behaved expansion in powers of k 2 only to order k 2 . Thus 
the scattering length and the effective range defined by Eq. (TTTT) are the only 
well-defined coefficients in the effective range expansion for atoms. 

The important energy scales for low-energy alkali atoms include the hyperfine 
splitting -Ehyperfine; the magnetic energy scale [iB, and the van der Waals energy 
scale defined by 



E vdw = (m'C/h 6 )- 1 / 2 . (18) 

For Rb atoms, the van der Waals energy scale is -E v dw = 6.4 x 10~ 9 eV. This 
is more than 3 orders of magnitude smaller than the hyperfine splitting E^. 
The magnetic energy scale fiB is comparable to -E v dw when the magnetic 
field is about 0.54 Gauss. The van der Waals energy scale corresponds to a 
temperature of 7.4 x 10 -5 K. We consider atoms to be ultracold if their kinetic 
energies are much smaller than the van der Waals energy scale. In subsequent 
sections, we will review scattering models that can describe such ultracold 
atoms. 
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2.4 Large scattering length 



The natural scale for the scattering length a of a specific hyperfine state is the 
van der Waals length defined by 



4dw = (mC 6 /ft 2 ) 1/4 . (19) 

A scattering length can be orders of magnitude larger than £ v< iw only if there 
is a fine-tuning of some parameter in the Hamiltonian. If \a\ is much larger 
than £ v dw, we refer to (X ctS 8b large scattering length. If a is large, ultracold 
atoms have an important energy scale that is much smaller than the van der 
Waals energy scale -E v dw, namely h 2 /(ma 2 ). 

For particles with short-range interactions, the scattering length a is consid- 
ered to be large if its absolute value is much larger than the range. Particles 
with large scattering lengths have universal low-energy properties that depend 
on the scattering length but are insensitive to other details of the interactions. 
The known universal properties have been summarized in a recent review [2]. 
In the case of atoms, the relevant range is the van der Waals length scale £ v dw 
defined in Eq. (1191) . In the 2-atom sector, the universal properties are rather 
simple. The cross section for elastic scattering of bosonic atoms in the same 
spin state with energy in the region \E\ <C E vdW , where -E v dw is the van der 
Waals energy scale defined in Eq. (TT81 . is given by the universal formula 



, v 8ira 2 . . 

aelastic(£) = — 2 = . (20) 

1 + ma z hi 



Shallow bound states whose binding energies are much smaller than -E v dw are 
also universal. If a < 0, there are no such bound states. If a > 0, there is a 
single shallow bound state with binding energy 



E D = l/(ma 2 ). (21) 
We will refer to this state as the shallow dimer. 

The natural scale for the scattering lengths a s and a t is the van der Waals 
length defined in Eq. (fT9l . The scattering length a s or a t can be much larger 
than £ vdw only if the depth of the potential V s (r) or V t {r) is tuned so that there 
is a bound state very close to the threshold. If a s or a t is large, several of the 
scattering lengths af mf of the hyperfine channels will also be large. Examples 
of alkali atoms that have large scattering lengths at zero magnetic field are 
6 Li (a t « -2160 a ), 85 Rb (a. » +2800 a ), and 133 Cs (a t » +2400 a ). 
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A large scattering length a/ m/ in a specific hyperfine channel can be obtained 
by tuning the magnetic field B to a Feshbach resonance Bi. If the magnetic 
field is near B iy there is a diatomic molecule near the scattering threshold for 
two atoms in the \f,mf] B) hyperfine state. The diatomic molecule is a bound 
state of atoms in a spin state with a higher threshold. Since the magnetic 
moment of the diatomic molecule differs from twice the magnetic moment 
of the \f, mf\ B) atom, the energy of the molecule can be made to cross the 
scattering threshold by changing the magnetic field B. The Feshbach resonance 
occurs at the value Bi for which the molecule is at resonance with two atoms 
at the scattering threshold. The existence of Feshbach resonances in Cs atoms 
was first pointed out by Tiesinga, Verhaar, and Stoof [3]. The first Feshbach 
resonances to be observed were the resonances at 853 G and 907 G for 23 Na 
atoms in the |1, —1) hyperfine level [I] and the resonance at 155 G for 85 Rb 
atoms in the |2, —2) hyperfine level [5]. The physics of Feshbach resonances in 
cold atoms has been summarized in a recent review [BJ. 

At a generic value Bq of the magnetic field, the scattering length a(B) for a 
specific hyperfine channel varies slowly as a function of the magnetic field B. 
It can therefore be expanded as a power series around B = B : 



At a Feshbach resonance, the scattering length is not a smooth function of 
B. When the magnetic field is well below the Feshbach resonance Bi, the 
scattering length has some off- resonant value a^ g . As B increases through the 
value Bi, a(B) diverges to +oo or — oo, jumps discontinuously to — oo or +oo, 
and then approaches the original off-resonant value db g . The scattering length 
must pass through zero at some value of B, either before or after Bi. If we 
denote that value by Bi + A i; the scattering length can be expressed as 



where a-^ g {B) is a smooth function of the magnetic field that can be expanded 
as a power series around Bi as in Eq. (|22|) . The parameters Bi and A, are 
convenient, because it is often possible to measure the locations of divergences 
and zeroes of the scattering length more accurately than its value a(B) at any 
particular value of B. 



a(B) = a(B ) + a'(B )(B - B ) + . . . . 



(22) 




(23) 
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3 Scattering models 



The fundamental theory described in Sec. [2] provides an extremely accurate 
description of systems consisting of atoms in any of the hyperfine states as long 
as their energies are small compared to the electronic excitation energy of the 
atom. If we restrict our attention to ultracold atoms, whose energies relative to 
a specific scattering threshold are small compared to the van der Waals energy 
scale -E v dw defined in Eq. (TT5I) . we can describe the system accurately using a 
much simpler model. The effects of virtual atoms whose energies are far from 
the scattering threshold can be taken into account through the parameters of 
the model. 

We take the atoms of primary interest to be those in a specific hyperfine state 
\f,rrif]B). We choose the zero of energy to be the scattering threshold for 
atoms in that hyperfine state. We restrict our attention to ultracold atoms 
whose energies satisfy \E\ <C -E v dw- If ^B <C -EVdw, this energy constraint can 
be satisfied by atoms in any of the 2/ + 1 spin states in the same hyperfine 
multiplet. For larger values of B, only atoms in the single hyperfine state 
\f,mf]B) satisfy this energy constraint. 



3.1 Potential models 

One class of models that can be used to describe the behavior of ultracold 
atoms in the spin state of interest is potential models. A potential model is 
specified by interaction potentials between the constituents of the model. In 
the simplest potential models, the only constituents are atoms in the spin 
state of interest. In this case, the potential model is specified by a single 
potential V(r). The potential V(r) must be tuned so that the solutions to 
the single-channel Schrodinger equation for that potential give the same re- 
sults for low-energy scattering cross sections and low-energy bound states as 
the solution to the coupled-channel Schrodinger equation for the fundamental 
theory. The potential V(r) need not resemble either of the potentials V s (r) 
or V t (r) in the fundamental theory as long as it accurately reproduces low- 
energy observables for atoms in the spin state of interest. If we restrict our 
attention to sufficiently low energies that the scattering length a is the only 
relevant interaction parameter, then all that is required is that the potential 
be short-ranged and give the correct scattering length. 

A potential model can be formulated as a quantum field theory. If the only 
constituents in the model are atoms in a single hyperfine state \ f, mf, B), there 
is a single quantum field i[)(r) which annihilates atoms in that hyperfine state. 
The Hamiltonian is 
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d 3 r Vip^ ■ Vip 



2m 




(24) 



The time- dependent quantum field ip{r,t) satisfies equal-time commutation 
relations. If the atom is a boson, the commutation relations are 



A more elaborate potential model could describe N—l additional spin states of 
the atoms. The model would be specified by the diagonal potentials for each 
pair of spin states and by coupling potentials between pairs of spin states. 
These potentials need not resemble those of the fundamental theory as long 
as the solutions to the Schrodinger equation for all the coupled channels re- 
produce the predictions of the fundamental theory for low-energy observables 
involving the spin state of interest. 



3.2 Scattering Models 

Another class of models that can be used to describe the low-energy behavior 
of the atoms is scattering models. A scattering model can be specified by 
the T-matrix elements for scattering of the constituents of the model along 
with any additional parameters that are required to make the model well- 
defined. In the simplest scattering models, the only constituents are atoms 
in a single spin state. At sufficiently low energies, the atoms have only S- 
wave scattering and the T-matrix element reduces to a function T(k) of the 
magnitude of the relative momentum k only. Since it is specified in terms of 
the low-energy observables T(k), a scattering model provides a particularly 
natural description for ultracold atoms. 

A more elaborate scattering model could include other spin states of the atom 
as additional degrees of freedom. In this case, the T-matrix elements for all the 
pairs of spin states could be part of the specification of the model. In contrast 
to potential models, diatomic molecules can also be included as explicit degrees 
of freedom in a scattering model. In this case, T-matrix elements involving the 
diatomic molecules could also be part of the specification of the model. 





ip{r,t),ip j (r',t) =5 3 (r-r'). 



(25b) 
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3.3 Hamiltonian quantum field theories 



A scattering model can be formulated as a quantum field theory. If the T- 
matrix element T{k) that is part of the specification of the scattering model is 
an analytic function of k 2 at k = 0, the T-matrix element can be reproduced 
by a local quantum field theory. A local quantum field theory is one whose 
interactions only involve the quantum fields and their derivatives at the same 
point in space. If a local quantum field theory has a Hamiltonian formulation, 
its Hamiltonian can be expressed as the integral over space of a Hamiltonian 
density: 



H = J d 3 r(H hcc + H mt ). (26) 

If the quantum field tp(r) that annihilates atoms in the spin state of interest is 
the only quantum field in the model, the free term in the Hamiltonian density 
is 



n ircc = ^-V^ -V^j. (27) 
2m 

The interaction term in the Hamiltonian density must be a polynomial in ip 
and ip^ and their gradients, such as 



+ i Po 



tV2 



W> • W> + 2^tyv^ • W> + i? V^ t • W> f ) + . . . . (28) 



One advantage of a Hamiltonian quantum field theory is that it guarantees 
that the space of quantum states is a Hilbert space. Since the Hamiltonian 
H is a hermitian operator, its eigenstates \n) form a complete set of states. 
The completeness relation can be expressed formally as J2n |n)(n| = 1. If we 
express the S-matrix in the form S — 1 + iT, the completeness relation can 
be used to express the unitarity condition S^S = 1 in the form 



2lm(i\T\i) = Y / \(n\T\i)\ 2 . (29) 

n 

This form of the unitarity condition is called the Optical Theorem. Note that 
the right side is the sum of positive-definite terms. 
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3.4 Lagrangian quantum field theories 



It is sometimes useful to consider a more general class of quantum field theories 
that have a Lagrangian formulation but do not necessarily have a Hamiltonian 
formulation. The Lagrangian for a local quantum field theory can be expressed 
as the integral over space of a Lagrangian density: 



L = Jd 3 r (£ frec + £ int ) . (30) 

If the quantum field ip{r) that annihilates atoms in the spin state of interest is 
the only quantum field in the model, the free term in the Lagrangian density 
is 



£ *-=5 i (*'s*-s* , *)-s; v * t - v *- (31) 

This form is manifestly real. If we ignore terms that are total time derivatives 
or total divergences, it can be written in a more compact form: 



£ - = (4 + 3 *■ < 32) 

If the quantum field theory has a Hamiltonian formulation, the interaction 
term in the Lagrangian density is simply £ mt = However, Lagrangian 

quantum field theories are a more general class of models than Hamiltonian 
quantum field theories. They can, for example, have interaction terms that 
involve time derivatives, such as 



C mt = ImX'^y (i| + ^ ^ + . . . . (33) 

More generally, £i n t must be a polynomial in ip, their gradients, and their 
time-derivatives. 

If all we have is a Lagrangian formulation of the quantum field theory, there 
is no guarantee that the Optical Theorem will have the form in Eq. (T29|) . 
where Im (i\T\i) is expressed as the sum of positive-definite terms. The general 
expression may have the form 



2Im(i|T|i) = ^|HT|2>| 2 -Om|T|2)| 2 , ( 34 ) 
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where the states labeled \n) and \m) together form a complete set. One way to 
interpret such a relation is that the space of quantum states is a complex vector 
space with indefinite norm. The completeness relation for this vector space has 
the form J2 n \ n )( n \ ~ E m l m )( m l = 1- The states \m) can be interpreted as 
states with negative norm. 

The problem of negative-norm states arises in covariant formulations of Quan- 
tum Electrodynamics. However, one can use gauge invariance to show that 
there is a physical subspace of the complex vector space spanned by the 
positive-norm states |n) such that if \i) is in the physical subspace, Im (i\T\i) 
can be written as the sum of positive-definite terms: 

2Im(i|T|i) = Y, K™l r l*)| 2 - ( 35 ) 

physical n 

This requires a cancelation between positive and negative terms in Eq. f)34p . 
In QED, this cancelation is guaranteed by gauge invariance. In a scattering 
model for atoms, there is no obvious mechanism to guarantee such a cancela- 
tion. Thus the existence of negative-definite terms in the imaginary part of a 
diagonal T-matrix element is more problematic. 



3. 5 Symmetries 

The symmetries of the fundamental theory can be used to constrain the Hamil- 
tonian or the Lagrangian of a local quantum field theory. We will give the 
constraints only for the simple case of a single quantum field ip(r,t): 

• phase symmetry. The phase transformation is tp(r,t) — > e i/)(r,t). Invari- 
ance under this transformation requires that every term in 7ii nt or £; nt have 
an equal number of factors of ip and ipK 

• translational symmetry. This is implemented automatically in a local quan- 
tum field theory for which Ti int or £ int is a function of ijj, ijj^, their gradients, 
and their time-derivatives. 

• time translational symmetry. This is also implemented automatically in a 
local quantum field theory. 

• rotational symmetry. If the model involves only a single quantum field 
ip(r,t), a rotation acts only on its vector argument r. The requirement 
of rotational symmetry is that 7ii nt or £ int be a scalar under rotations. This 
requires each gradient V* to have its index i contracted with that of another 
gradient V* or with a Levi-Civita tensor e ljfc . 

• parity symmetry. This can be used to exclude terms in 7i int or £ int with a 
Levi-Civita tensor. 
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Galilean symmetry imposes strong constraints on the interaction terms in H m t 
or £; nt . The Galilean transformation with velocity vector v in the fundamental 
theory is given in Eq. (1121) . The corresponding Galilean transformation for the 
field ip{r,t) is 

ip(r, t) — ► e ^r-i(^ 2 )t^^ r _ vty t y ( 36 ) 

One way to construct interaction terms that are Galilean invariant is to con- 
struct them out of building blocks with simple behavior under Galilean trans- 
formations. One simple building block is the product which is Galilean in- 
variant. Thus terms constructed out of ip^ip and its gradients, such as V 1 (i/j^i/j) 
and V'V^ty), 

are Galilean invariant. The first term in 7i int in Eq. (I28p is 
Galilean invariant, because it can be expressed as (ip^ip) 2 . The second term 
in Tii-at is a ls° Galilean invariant, because it can be expressed in the form 

v(v>ty) ■ v(v>ty). 

Another simple building block is ^Vt, where the operator V is defined by 

*VV = X(V 2 ^) - 2(V X ) • (V^) + (V 2 x)^. (37) 

The operator ipS/ 2 ^ transforms under Galilean transformations in the same 
way as ip 2 . A further set of simple building blocks is the combination [ij^ + 
V 2 / ' (2Nm)}ip N , where iV = 1, 2, . . .. They transform under Galilean transfor- 
mations in the same way as ij) . This combination with N = 1 appears in the 
expression for £f ree in Eq. ( 1321) . This combination with N = 2 appears in the 
interaction term in Eq. fl33l) . 

If the energy restriction allows atoms in multiple spin states related by a 
symmetry, the symmetry can be imposed on the corresponding quantum fields 
VVnty, t). For example, if B = 0, the 2/ + 1 spin states in a hyperfine multiplet 
form an irreducible (2/+l)-dimensional representation of the SU(2) symmetry 
group of angular momentum. There is also a f/(l) symmetry corresponding 
to multiplying all the fields if) m (r,t) by the same phase. 

3.6 Two-body observables 

In a quantum field theory, physical observables can be conveniently encoded in 
the correlation functions of the theory. The physical observables in the 2-atom 
sector can be encoded in the amputated connected Green's function A for two 
atoms in the asymptotic past to evolve into two atoms in the asymptotic 
future. In general, this amplitude is a function of the energies E\ and E 2 and 
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the momenta pi and P2 of the two incoming atoms and the energies E[ and 
E' 2 and the momenta p[ and p' 2 of the two outgoing atoms. The atoms can 
be off their energy shells: Ei need not be equal to p 2 /{2m). The total energy 
E and the total momentum P are conserved: E 1 + E 2 = E[ + E' 2 = E and 
Pi + P2 = Pi + p'i = P- Galilean invariance constrains the amplitude to be 
a function of E — P 2 /(4m) and the relative momenta k = \{pi — P2) and 
k' = \{p'i — p 2 ). We can simplify the problem by working in the center-of- 
mass frame P = 0. If the energy is restricted to be low enough that there 
is S-wave scattering only, the amplitude in the center-of-mass frame depends 
only on the magnitudes of k and k', so it reduces to a function A(E; k, k') of 
three variables. 

The amplitude A(E; k, k') encodes all information about low-energy 2-atom 
scattering. The T-matrix element T(k) for the S-wave scattering of two atoms 
with relative momentum k into two atoms with relative momentum k' is 
obtained by evaluating A(E; k, k') at the on-shell point obtained by setting all 
the atoms on their energy shells: E = 2{k 2 /2m) = 2(k' 2 /2m), which requires 
k — k'. The T-matrix element is 



The cross section for the elastic scattering of two identical bosonic atoms with 
total kinetic energy E = k 2 /m is obtained by squaring the T-matrix element 
and multiplying by the phase space factor mk/iAir) and by the flux factor 
m/{2k): 



Information about S-wave bound states is also encoded in A(E; k, k') through 
its poles in the energy E. The amplitude A(E; k, k') is a double- valued function 
of the energy variable E, which takes its values on the two-sheeted complex 
plane with a branch point at E — and a branch cut along the positive E 
axis. The interpretation of the pole in E depends on where the pole resides 
in the complex E plane. For simplicity, we assume that the atom has no spin 
states with lower energy. A pole on the negative real axis of the physical sheet 
of the complex energy E corresponds to a stable bound state. A pole on the 
negative real axis of the unphysical sheet is called a virtual state. A pole on the 
unphysical sheet just below the positive real axis corresponds to a scattering 
resonance. 



T{k) = A(k 2 /m;k,k). 



(38) 



a(E) 




\T(k)\ 2 . 



(39) 
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3. 7 Renormalization 

One complication of local quantum field theories is that they are inherently 
singular at short distances. The theory is well-defined only if it is modified in 
some way to suppress the effects of particles that approach arbitrarily close to 
each other or that have arbitrarily high energy. A convenient way to do this 
is to impose an upper limit on the wavevector of virtual particles: \k\ < A. 
The variable A is called the ultraviolet cutoff. The singularities then appear 
in the limit A — > oo. However, it is possible to avoid the singularities by using 
a procedure called renormalization. The parameters of the model, such as A 
and po m Eq. (|28|) . are made functions of A in such a way that the model has 
a well-behaved and nontrivial limit as A — > oo. 

There are two rather different approaches to renormalization. The simplest 
approach conceptually is called G ell- Mann- Low renormalization. This is the 
approach that was followed in the original development of renormalized pertur- 
bation theory for Quantum Electrodynamics. A quantum field theory model 
can be specified by a Hamiltonian density with a finite number of interaction 
terms: 

N 

H = H {ree + ]T c n (A)O n . (40) 

n=l 

If the coupling constants c„(A) can be tuned as functions of the ultraviolet 
cutoff so that physical observables have well-behaved and nontrivial limits 
as A — > oo, the model is called renormalizable. It may be necessary to add 
operators to the set {Ox, . . . , On} or remove some of the operators to make 
the model renormalizable. In nonrelativistic quantum field theories in which 
the total number of particles cannot change, the physical observables can be 
classified according to the particle number. If the coupling constants can be 
tuned so that all A-particle observables have well-behaved limits as A — > 
oo, the model is called renormalizable in the N-atom sector. If the model is 
renormalizable in the A-atom sector, it is renormalizable in the n-atom sector 
for all n < N . In this paper, we will for the most part consider only the 
constraints from renormalizability in the 2-atom sector. 

The other approach to renormalization is called Wilsonian renormalization. 
This approach was followed by Ken Wilson when he developed the renormal- 
ization group to understand strongly-interacting quantum field theories. One 
begins by considering the most general model consistent with the symme- 
tries. The Hamiltonian density in Eq. (piOl is generalized to one with infinitely 
many operators O n , each with its own coupling constant c„(A). The free term 
?if re e is included as just another operator with a coupling constant Cf ree (A). 
We define the engineering dimension d n of an operator O n by specifying the 
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dimensions of if) and V to be | and 1, respectively. For this purpose, the mass 
m is treated as a dimensionless conversion factor. The engineering dimension 
of the free term in Eq. (1271) is then 5. An A-atom operator with N factors of 
if), N factors of if)', and D factors of the gradient V has engineering dimension 
3 A + D. It is convenient to introduce dimensionless coupling constants c n (A) 
defined by 



c„(A) = A d "- 5 c n (A). (41) 

The condition that low-energy observables defined by the Hamiltonian be in- 
dependent of A defines a flow on the infinite-dimensional space of dimension- 
less coupling constants c n . This flow can have a fixed point whose coordinates 
c n are independent of A. The corresponding fixed-point Hamiltonian Tt* de- 
scribes a scale- invariant system. Near the fixed point, a convenient basis of 
operators consists of scaling operators 0* n that have the property that when 
their coupling constants are infinitesimally small, they scale as definite powers 
of A: 



Cn(A)~A^. (42) 

We refer to the exponent 7„ as the scaling dimension of the operator 0* n . An 
operator with scaling dimension j n < is called relevant, because it becomes 
increasingly important as A decreases towards and thus has a dramatic 
effect on the physics at low energies. An operator with scaling dimension 
7„ > is called irrelevant, because it decreases in importance as A decreases 
towards 0. An operator with scaling dimension j n = is called marginal. The 
Hamiltonian density can be expanded around the fixed point in terms of the 
scaling operators: 



H = H* + £ c n (A)0* n . (43) 

n=l 

Once the scaling operators have been identified and their scaling dimensions 
7„ determined, the model can be truncated to include only a finite number 
of terms in the Hamiltonian. A consistent truncation must include all the 
operators that are relevant and marginal. The behavior of the system near the 
fixed point can be described more accurately by including irrelevant operators. 
The accuracy of the description can be systematically improved by adding 
operators with increasing values of j n . 

The simplest example of a renormalization group fixed point is the free theory 
with the Hamiltonian density 7if ree in Eq. (I2TI) . This fixed point is called the 
trivial fixed point. The forms of the scaling operators depend on the choice 
of ultraviolet cutoff or regularization. With a scale-invariant regularization 
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scheme, such as dimensional regularization, the scaling operators consist of 
terms that all have the same engineering dimensions. A scaling operator with 
N factors of tp, N factors of if)' , and D factors of V has scaling dimension 
3iV + D — 5. The 2-body interaction terms in Eq. (I28p are irrelevant operators 
with scaling dimensions 1 and 3, respectively. The simplest 3-body interaction 
term (^) 3 -?/> 3 is an irrelevant operator with scaling dimension 4. The successive 
truncation of the Hamiltonian to include operators with increasingly higher 
scaling dimension corresponds to a combination of the gradient expansion and 
an expansion in the number of particles. 

The system of nonrelativistic particles with short-range interactions also has a 
nontrivial fixed point. In terms of the effective range expansion in Eq. (IT7|) . the 
fixed point corresponds to the limit in which the scattering length a is infinite 
and the effective range r s and all higher coefficients in the expansion are 0. The 
dimensions of the scaling operators at the fixed point were deduced by Kaplan, 
Savage and Wise [7115] and by van Kolck [5] . There is a relevant operator with 
scaling dimension —1 which if added to the fixed-point Hamiltonian gives a 
model with a finite scattering length. There are also irrelevant operators with 
scaling dimensions 1, 3, 5, .... 



4 Zero-Range Model 

The simplest scattering model for ultracold atoms is the Zero-Range Model, 
whose only parameter is the scattering length a. The Zero-Range Model can 
be defined by the truncation of the effective range expansion in Eq. (fTTI) after 
the constant term: 

kcotd {k) = -1/a. (44) 

The effective range and all the higher coefficients in the effective range expan- 
sion in Eq. (ITTj) are exactly 0. 

The Zero-Range Model provides a natural description for atoms with a large 
scattering length that arises from the accidental fine tuning of the depth of the 
interatomic potential. In such the scattering length typically changes 

very slowly with the magnetic field B. Its dependence on B can be taken 
into account by allowing the parameter a in Eq. (|44p to be a smooth function 
of B. The Zero- Range Model also provides a minimal description for atoms 
with a large scattering length that arises from any other mechanism. If that 
mechanism is a Feshbach resonance, the dependence of the scattering length 
a on B has the form given in Eq. (j23p . 
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4-1 Hamiltonian 



The Zero-Range Model can be described by a quantum field theory with the 
single complex field ip. The Hamiltonian density is the sum of the free term 
in Eq. (f2~T|) and the interaction term 



ttint = ^Ao(VW- (45) 

The interaction term corresponds to a momentum-independent vertex with 
the Feynman rule —2 An. Such an interaction is singular, so the model requires 
an ultraviolet cutoff A on the momenta of virtual atoms. 

The Zero-Range Model is renormalizable in the 2-body sector. As shown in 
Sec. I4.2[ the bare coupling constant Ao in Eq. (j4"5l) can be adjusted as a function 
of A so that 2-body observables do not depend on A. The scattering length 
has the desired value a if the bare coupling constant is 



VAHfl-V)" 1 ^. (46) 

This expression can be inverted to express the scattering length as a function 
of the bare coupling constant and the cutoff: 



m 
8^ 



1 m 



— + — A 

o 



An 4tt 2 ' 



-i 



(47) 



The expression for the bare coupling constant in Eq. (l4*6|) provides a simple 
illustration of Wilsonian renormalization, which was discussed in Sec. 13.71 The 
operator (ip^ip) 2 in Eq. (1451) has engineering dimension 6. The corresponding 
dimensionless coupling constant is therefore AAn(A). This dimensionless cou- 
pling constant has fixed points for two values of the scattering length: a = 
and a = ±oo. The trivial fixed point a = corresponds to a noninteracting 
theory. For infinitesimal values of a, the dimensionless coupling constant is 
AAn(A) « 8iraA/m. Thus near the trivial fixed point, (ip^ip) 2 is an irrelevant 
operator with scaling dimension 1. The nontrivial fixed point a = ±oo is com- 
monly referred to as the unitary limit. At the fixed point, the dimensionless 
coupling constant is Aq ~ — 47r 2 /m. For infinitesimal values of 1/a, the de- 
viation of the dimensionless coupling constant from its fixed-point value is 
AAn(A) — Aq ~ — 27r 3 /(maA). Thus near the unitary fixed point, (ip^ip) 2 is a 
relevant operator with scaling dimension —1. 
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X-XX'XxX- 

Fig. 3. The series of diagrams whose sum is the amputated connected Green's 
function for atom-atom scattering in a model with contact interactions only, such 
as the Zero-Range Model. 

4-2 Green's function 




The amputated connected Green's function for 2 atoms to evolve into 2 atoms 
can be calculated analytically by summing the series of diagrams shown in 
Fig. [31 Each vertex gives a factor of — i\ . Each loop gives the product of two 
propagators, which must be integrated over the momenta and energies of the 
virtual atoms. In the center-of-mass frame, the integral is 



du f d 3 k % i 

■ 10 — k 2 1 (2m) + ie E — u — k 2 / (2m) + ie 

(4£ 



(2tt) 3 E -k 2 /m + ie' 

where E is the total energy in that frame. The energy integral has been evalu- 
ated using contour integration. The remaining momentum integral is ultravio- 
let divergent. It can be regularized by imposing a momentum cutoff: |fc| < A. 
For each loop, there is also a symmetry factor of |. The resulting expression 
for the amplitude is 



iA{E) = -i\ £ 

n=0 

The amplitude depends on the total energy E of the two atoms in the center- 
of-mass frame, but not on the momenta of the incoming or outgoing atoms. 
The sum of the geometric series in Eq. (|49|) is 



H'Aq 



1 



d 3 k 



2 J (2tt) 3 E — k 2 /m + ie 



(49) 



A(E) = -Ac ^1 - fI (E)j , (50) 

where the integral Iq(E) is defined in Eq. ( 1A.1I) . The expression for the integral 
is particularly simple if we take A to be so much larger than {m\E\) 1 l 2 that we 
can neglect terms that decrease as inverse powers of A/ \m\E\) l l 2 as A — > oo. 
The integral is given in the Appendix in Eq. (1A.3I) . The resulting expression 
for the Green's function is 
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Fig. 4. The integral equation for the amputated connected Green's function for 
atom-atom scattering in a model with contact interactions only, such as the Ze- 
ro-Range Model. 



A(E) 



A + 4vr 2 



m . m 

A K 

8tt . 



(51) 



where k is proportional to the square root of the energy: 



f —mE — ie 



,1/2 



(52) 



The amplitude in Eq. (15ip can also be obtained by solving the Lippmann- 
Schwinger integral equation, which is represented diagrammatically in Fig. H] 
and can be written 



iA(E) = -i\ Q + (-iX )-J 



d 3 k 



(2vr) 3 E - k 2 /m + ie 



(iA{E)). 



(53) 



Since A(E) does not depend on the momenta of the incoming or outgoing 
particles, this is just a linear equation for A(E) and its solution is Eq. (l5~Tj) . 

The dependence of the amplitude in Eq. ( |5"TI) on the ultraviolet cutoff can be 
eliminated by choosing the bare coupling constant Ao(A) to have the form in 
Eq. f)46p . where a is an arbitrary parameter with dimensions of length. The 
renormalized expression for the amplitude is then obtained by using Eq. fPfSl 
to eliminate Aq in favor of a: 
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A(E) = 



m 1 — an 



(54) 



Since A(E) encodes all the physical 2-body observables, we conclude that the 
Zero-Range Model is renormalizable in the 2-body sector. 



A simple way to reproduce the renormalized expression for A(E) in Eq. ([5] 
is to use dimensional regularization, which automatically subtracts power ul- 
traviolet divergences such as the linear term in A in Eq. flA.3h . With this 
regularization scheme, the integral Io{E) is given simply by Eq. flA.6aj) . Since 
no renormalization is necessary in this case, we call the parameter A instead 
of Aq. The amplitude in Eq. (150]) then reduces to 



A(E) 



' 1 m 

.A ~ 8^ K 



n -1 



(55) 
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This agrees with Eq. if we set A = 87ra/m. 



4-3 T -matrix dement 

The T-matrix element for atom-atom scattering with relative momentum k 
can be obtained from A(E) in Eq. (I54j) by setting E equal to the total kinetic 
energy 2(k 2 /2m) of the two atoms: 

T(k) = -^—^-r. (56) 
m 1 + iafc 

By taking the low-energy limit k — > and comparing with Eq. (fl4l) . we identify 
a as the scattering length. Comparing with Eq. (fT6l) . we obtain the simple 
equation for the phase shift in Eq. (PHI) . 



^.^ Bound state 

The amplitude -4(-E) in Eq. fl54|) has a pole in the energy variable k defined 
in Eq. fl52|) . The pole is at k = 1/a. If a > 0, this corresponds to a pole in 
the energy at E = — 1/ (ma 2 ), which implies that there is a stable bound state 
with binding energy l/(ma 2 ). This bound state is the shallow dimer, whose 
binding energy is given in Eq. (T2"TT) . If a < 0, the pole at k = 1/a corresponds 
to pole at E = — 1/ (ma 2 ) on the second sheet of the complex variable E. Such 
a pole corresponds to a virtual state. 



4-5 Optical Theorem 

If the energy E is real, the imaginary part of the amplitude in Eq. (I5il) is 

. , i (87i/m)a 2 k „. „ s 167r 2 /m 2 .. 9 .. „. . , , 

Im^(E) = \ 7 ; 19 0(g) + '- 5{E + 1 / ma 2 61(a). 57 

1 + a^/r a 

The first term is nonzero for positive values of E = k 2 /m. The second term, 
which contributes only if a > 0, is a delta function at the negative energy 
—Ed- Note that the expression for ImA(E) in Eq. (IBTj) is positive definite, in 
accord with the Optical Theorem in Eq. (|29j) . This implies that there are no 
negative-norm states in the 2-atom sector of the Zero- Range Model. 
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4-6 Three-body sector 



In the 2-body sector, the Zero-Range Model is completely defined by a sin- 
gle parameter: the scattering length a. In the 3-body and higher sectors, this 
condition is not sufficient to specify the model completely because of the Efi- 
mov effect [TOlfTTj . In the unitary limit a — > ±00, there are infinitely-many 

(n) 

arbitrarily-shallow trimers called Efimov states whose energies — have an 
accumulation point at the 3-atom threshold. For a = ±00, the spectrum is 



R 



(n) 



m 



(5* 



where n # is an integer, has dimensions of momentum, and Sq is the tran- 
scendental number that satisfies 



s cosh — = — sinh — . (59) 

Its numerical value is So ~ 1.00624. The spectrum in Eq. (I58p defines a 3- 
body parameter /«* that is independent of a. One can interpret Ti 2 Kl/m as 
the approximate binding energy of the trimer labeled n*. A different choice 
of gives a value of k* that differs by an integer power of e^ ~ 22.7. 
The ratio of successive binding energies in Eq. (158]) is e 27r//s ° ~ 515.03. This 
spectrum in Eq. (|58|) is consistent with an asymptotic discrete scaling sym- 
metry with discrete scaling factor ~ 22.7. The 3-body observables are 
uniquely determined only if one specifies both the scattering length a and the 
3-body parameter k*. The 3-body parameter be adjusted by adding a 

(ip^i/j) 3 term to the Hamiltonian [T2"][T3"] . The renormalization of its coefficient 
is governed by a limit cycle instead of a fixed point [2]. Thus if the Zero-Range 
Model is extended to include a (ip^ijj) 3 interaction term, it is renormalizable 
in the 3-atom sector. 

It is an open question whether the Zero-Range Model is renormalizable in 
the 4-atom and higher sectors. One problem is that the spectrum of Efimov 
states in Eq. fl58|) is unbounded from below in the limit A — > 00, where A 
is the ultraviolet cutoff. Thus there can be no stable 4-atom bound states in 
the Zero-Range Model, because they can decay into an Efimov state and an 
atom. In spite of this pathology of the model, it is still possible that low-energy 
observables in the 4-atom and higher sectors have well-defined universal limits 
as A — ► 00 that depend only on the parameters a and k*. There is numerical 
evidence that this is the case in the 4-atom sector [T4]fT5] . 
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5 Effective Range Model 



The simplest generalization of the Zero-Range Model is the Effective Range 
Model, whose two parameters are the scattering length a and the effective 
range r s . The Effective Range Model can be defined by the truncation of the 
effective range expansion in Eq. (fTTl) after the k 2 term: 



kcot5 {k) = -l/a + \r s k 2 . (60) 

All the higher coefficients in the effective range expansion are 0. The Effective 
Range Model provides a natural description for atoms that have a large scat- 
tering length and a large effective range because of a double fine-tuning. It 
can also describe atoms with a large scattering length and a natural effective 
range more accurately than the Zero-Range Model. 

Phillips, Beane, and Cohen showed that the phase shift in Eq. (|6"U|) can be 
reproduced by a renormalizable local Hamiltonian quantum field theory with 
a single atom field only if the effective range r s is negative [16] . This will be 
verified below in Sec. 15.31 In Sec. [8], we will show that the phase shift in Eq. (1601) 
can also be reproduced by a renormalizable local Hamiltonian quantum field 
theory with an atom field and a molecular field, but again only if r s < 0. In 
Sec. 16. 3[ we will show that the phase shift in Eq. (1601) with either sign of r s can 
be reproduced by a Lagrangian quantum field theory with a time- derivative 
interaction. However, if r s > 0, the price that must be paid is a space of 
quantum states that includes negative-norm states. 



5. 1 Hamiltonian 



The simplest Hamiltonian quantum field theory that gives a nonzero effective 
range has the two interaction terms in Eq. (1281) . which can also be expressed 
in the form 



n int = \\^H) 2 + ^poV(^) ■ v^ty). (6i) 

The second interaction term can be expressed in an alternative form using 
Eq. ([37D: 



(62) 



Hint = ^A (^ 2 )V 2 - lp (^ 2 ) W</0 + (^)V 2 

4 16 
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The difference between the Hamiltonian densities in Eqs. ( IBTj) and (|62|) is a 
total divergence, so it gives zero when integrated over all space to get the 
Hamiltonian. The elementary vertex for particles with momenta pi and p 2 to 
scatter into momenta p[ and p' 2 is 



/' ( A + ^Po [(Pi - P2) 2 + (Pi - P2) 2 



Ao + po^ + A:' 2 ) , (63) 



where k = \\pi — P2I and fc' = — p' 2 | are the relative momenta in the 
initial and final states. 

As shown by Phillips, Beane, and Cohen, the Effective Range Model is renor- 
malizable in the 2-body sector provided the effective range r s is negative |16j . 
The bare coupling constants A and p i n Eq. (|62|) can be adjusted as func- 
tions of the ultraviolet cutoff A so that 2-body Green's functions and 2-body 
observables have finite limits as A — > 00. The intricate tuning that is required 
will be derived in Sec. 15.31 



Po(A) 
Ao(A) 



12tt 2 



mA 3 
m 



1 ± 



-12 ^ 



7rr Q A 



7T 



6 



Po(A) 2 A 5 + 



48vr 
mr s A 2 



2aA 

7T 



+ 



7rr s A 
12 



2aA nr<A 



+ ... 



(64a) 
(64b) 



Note that the dimensionless coupling constants AAo and A 3 p approach finite 
limits as A — >• 00. 



5.2 Green's function 



The amputated connected Green's function for two atoms can be calculated 
analytically by solving the Lippmann-Schwinger integral equation represented 
diagramatically in Fig. HI This equation can be solved most easily by assuming 
that the amplitude can be expressed in the form 



A(E; k, k') = A (E) + Ai(E) (k 2 + k' 2 ) + A 2 (E)k 2 k' 2 . (65) 

The integral equation then reduces to a system of linear equations for the three 
functions Aq(E), Ai(E), and A 2 (E) [IB] . The final result for the amplitude 
can be expressed in the form 



A{E; k, k') 



N(E;k, k') 
D{E) 



(66) 
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The numerator is a function of E, k, k', and A: 



N{E;k,k') = (\ 



m 
207 2 ' 



p^A 5 ) - 2p I 1 



m 
24^ 



-Po 
in 



1 + Y^PoA 3 ) (2k 2 + k 2 + k' 2 



i ^(A-|«)(« 2 + ^ 2 )(- 2 + /' i 



(67) 



where n = (—mE — is) 1 / 2 . The denominator is a function of E and A: 



D(E) 



m 
12^ 



■p A ; 



+ 



m 
4^ 



An 



m 
20^ 



-2p 1 



m 
24^2 



PoA d re 



A 



7T 



k \ . (6S 



Since the numerator is a linear combination of four independent functions of 
k, k, and fc', the dependence of the amplitude A(E] k, k') on A cannot be 
completely eliminated simply by taking the two bare parameters Ao and po to 
be functions of A. However, as shown by Phillips, Beane, and Cohen [16], if 
the bare parameters A and p are tuned as in Eqs. (1641) . the ratio in Eq. (1661) 
will have a finite limit as A — > oo. The renormalized amplitude obtained by 
this limiting procedure has the form 



A(E) = - 



8tt 



m 1 + ^r s aK 2 



CLK 



(69) 



It depends on the total energy E but not on the momenta k and 



5.3 T-matrix element 



The T-matrix element for atom-atom scattering with relative momentum k 
can be obtained from A(E) in Eq. (1691) by setting equal to the total energy 
k 2 /m of the two atoms: 



ml — ^r s ak z + za/c 



Comparing with Eq. (TT6I) . we obtain the equation for the phase shift in Eq. ([60 
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To drive the expressions for the tuning of the bare parameters in Eqs. (164|) . we 
start from the unrenormalized expression for the T-matrix element obtained 
by setting k! = k and k = —ik in the Green's function A(E; k, k') given by 
Eqs. flSSP, dSZD, and (ESJ): 



T(k) 



[1 + (m/12vr 2 )p A 
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[A - (m/20vr 2 )p 2 A 5 ] + 2 p [l + (m/24vr 2 )p A 3 ]£; 2 



+ 4^ A + ! ^j ■ (71) 

Expanding T(k)~ l in powers of and comparing with the effective range 
expansion in Eq. ( 1TT1) . we find that the scattering length and the effective 
range are 



m 



8vr [1 + (m/127r 2 )p A 3 ] 2 /[A - (m/20vr 2 )p 2 A 5 ] + (m/47r 2 )A ! 
32vr po[l + (m/127T 2 )p A 3 ] 2 [ 1 + (m/247r 2 )p A 3 ] 



ni 



[A - (m/207r 2 )p 2 A 5 ] ; 



(72a) 
(72b) 



Using Eq. (j72aj) to eliminate Ao in favor of a, the expression for the effective 
range in Eq. f!72bj) becomes 



_ 2mp [l + (m/247r 2 )p A 3 ] 
Ts ~ ~^ [1 + (m/127r 2 )p A 3 ] 2 

One can show that for the right side of Eq. ( 1731) to have a finite nonzero limit 
as A — > oo, the quantity 1 + (m/127r 2 )poA 3 in the denominator must approach 
OasA" 1 / 2 : 



A-^ 
2a 



(73) 



[1 + (m/12vr 2 )p A 3 ] 2 — (l - ^ + + . . ). (74) 

This equation implies that the bare parameter p must have the limiting be- 
havior given in Eq. f)64al) . 

If we insert the limiting expression in Eq. f!74|) into the expression for the 
scattering length in Eq. ( 172 aft , it reduces to 



m 
8tt 



-12/(7rr s A) 



A - (m/207r 2 )p 2 A 5 



7T 12 

1 "-T + T 

aA 7rr s A 



+ ... + 
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47T 5 
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(75) 
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This shows that the quantity Aq — (m/207r 2 )pQA 5 must approach as A 2 : 



m 2*5 4871 " ft 71 12 \ /t-rn\ 

A »-20^» A ^^( 1 -^A + ^ + -)' (76) 

It also implies that the bare parameter Ao must have the limiting behavior 
given m Eq. flBlbl) . 

After inserting the limiting expressions in Eqs. f)74p and (1761) . the T-matrix 
element in Eq. ( ITTj) reduces to 



l (k) = hm 

m a^oo 



-A ( 1 + — — - —j-k ) + -A + ik 

2aA 4A / 7T 



TT 



(77) 



Taking the limit A —>■ oo, we recover the renormalized expression for the T- 
matrix element in Eq. flTUl) . 

Since the left side of Eq. (1741) is positive definite, this equation implies that 
the effective range r s in this model must be negative. This result was first 
derived by Phillips, Beane, and Cohen [16J. They showed that this conclusion 
is more general. The effective range must be negative for any model with 
a Hamiltonian that includes gradient interactions of arbitrarily high order. 
This result is consistent with a proof that the self-adjoint extension of the 
Hamiltonian with a delta-function potential to a 2-parameter Hamiltonian 
with scattering length a and effective range r s can be made on a Hilbert space 
with positive definite norm only if r s < [17J. 



5.4 Bound states 

The amplitude A(E) in Eq. fl69|) has two poles in the energy variable k defined 
by Eq. fl52l) . The values of k at the poles satisfy a quadratic equation: 

1 - an + \r s aK 2 = 0. (78) 

If there is a positive real root Kj, there is a stable bound state with binding 
energy K 2 /m. If there is a negative real root Ki, the pole corresponds to a 
virtual state. The quadratic polynomial in Eq. (1781) has no real roots if a/r s < 
2. If a/r s > 2, it has two real roots: 



« ± = -(l±V 1 - 2r -/o)- (79) 
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The condition a/r s > 2 requires a and r s to either have opposite signs or to 
have the same signs and satisfy \a\ > 2\r s \. The root k_ is positive if a > 
and negative if a < 0. The root n + is positive if r s > and negative if r s < 0. 
Thus there can be 0, 1, or 2 stable bound states depending on the signs of a 
and r s . 

If \a\ is large compared to \r s \, the approximate solutions are k_ « 1/a and 
k + w 2/r a . The solution k_ corresponds to the shallow dimer if a > and to a 
shallow virtual state if a < 0. If r s < 0, as required by the renormalizability of 
the Effective Range Model, the solution k + corresponds to a virtual state. If 
r s > was allowed, the solution k + would correspond to a stable deeply-bound 
state. 



5.5 Optical theorem 

If the energy E is real and a/r s > 2, the imaginary part of the amplitude in 
Eq. flSHD is 



K+ 8{E + K 2 + /m)6{r s ) . 

(80) 

The first term is nonzero for positive values of E = k 2 /m. This is the only 
term if a/r s < 2. If a/r s > 2, there can also be delta function contributions 
at negative values of E. There is a delta function term with a positive coef- 
ficient if a > and another delta function term with a negative coefficient 
if r s > 0. Thus if r s > 0, the Optical Theorem includes a negative-definite 
term as in Eq. ( 134"|) . The bound state associated with this term can be inter- 
preted as a state with negative norm. The condition r s < for the absence of 
negative-norm states coincides with the condition for the renormalizability of 
the Effective Range Model. If there is a large scattering length \a\ ^> r s , the 
negative-norm state corresponds to a deeply-bound state. 



6 Other Single-Channel models 

In this sector, we describe briefly other scattering models that can be formu- 
lated as local quantum field theories with a single quantum field if}. 



Im A(E) 
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6.1 Higher- gradient interactions 



One straightforward generalization of the Effective Range Model is to include 
higher gradient terms in the interaction Hamiltonian. There are two indepen- 
dent terms with four gradients that give only 5*-wave scattering: 



The model whose Hamiltonian is the sum of 7if re e in Eq. (j27|) . TC- m t in Eq. flBTj) . 
and AHi n t in Eq. (IHTj) has four bare parameters that can depend on the ul- 
traviolet cutoff. If this model is renormalizable, one might expect it to be able 
to reproduce additional terms in the effective range expansion of the S"-wave 
phase shift in Eq. ( 1TT1) beyond the two terms that are shown. Such a model 
is not of much practical use for atoms, because the 1/r 6 tail of the van der 
Waals potential implies that fccot 5o(k) has an expansion in powers of k 2 only 
to order k 2 . 

In Ref. [18] . the Green's function A(E; k, k') for this model with an ultraviolet 
cutoff A was obtained by solving the Lippmann-Schwinger integral equation 
shown diagrammatically in Fig. |H The authors showed that the dependence of 
the T-matrix element on the ultraviolet cutoff could not be removed simply by 
choosing the bare parameters to be functions of A, as in the Zero- Range Model 
in Sec. HJ However, they did not investigate whether the model is renormal- 
izable in the sense of having a nontrivial limit as A — > oo, as in the Effective 
Range Model in Sec. Since the Effective Range Model is the special case 
/Uo = /^o = 0, it is clear that the four bare parameters can be chosen as func- 
tions of A so that the phase shift in the limit A — > 00 is Eq. ( 160]) with r s < 0. 
It is not known whether a more general phase shift can be obtained by taking 
such a limit. 

6.2 Dimensional regularization 

If there is a consistent way to renormalize a quantum field theory, one can often 
obtain the final result more simply by using dimensional regularization. In 
this regularization method, ultraviolet-divergent integrals over a 3-dimensional 
momentum are generalized to a variable number of dimension d, evaluated in 
a range of d for which they converge, and then analytically continued to d = 3. 
This regularization prescription treats power divergences very differently from 
logarithmic divergences. Power- divergent integrals are those that diverge as 




+// (^ 2 ) f (Vv 2 vV) + (Vv 2 v 2 ^) 



(81) 
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a positive power A p with a momentum cutoff A. Logarithmically-divergent 
integrals are those that diverge as log A as A — > oo. Dimensional regularization 
sets power divergence to 0, while logarithmic divergences appear as poles in 
d-3. 

In the Effective Range Model, which is defined by the interaction Hamiltonian 
in Eq. (IBTI) . the only ultraviolet divergences that appear are power divergences 
A p with p = 1,3, 5. Dimensional regularization sets these divergences to 0. 
Thus the Green's function A(E; k, k') with dimensional regularization can be 
obtained from the results in Eqs. (JHSD, (IBTI) . and (I6"%l) with a momentum cutoff 
simply by setting A = 0. The resulting expression is 



A(IPt us A + p(k 2 + k' 2 ) + (m/frr) A(* 2 + fc 2 )(^ 2 + k' 2 ) 
A V> k > k >- l-(m/87r)(A-2p«2)« • (8 } 

Since no renormalization is required, we have replaced the bare parameters Ao 
and po by A and p. The T-matrix element T(k) for atom-atom scattering can 
be obtained from the amplitude A(E; k, k') in Eq. (152]) by setting k! = k and 
K = —ik: 



T(/C) l/(A + 2pP)+2(m/8vr)A:' (83) 

Expanding T(k)^ 1 in powers of k and comparing with the effective range 
expansion in Eq. flTTJ, we find that the parameters A and p are related in 
a simple way to the scattering length and the effective range: A = Sna/m, 
p = 27ta 2 r s /m. If A and p are eliminated in favor of a and r s , the T-matrix 
element in Eq. ( 1831) reduces to 



This result corresponds to the S"-wave phase shift 



kcot5 (k) = j - T . (85) 

a 1 1 + ^r s ak z 



The phase shift in Eq. ( 1851) obtained using dimensional regularization has a 
different dependence on k from the phase shift in Eq. (|60|) obtained using 
a momentum cutoff regularization. Renormalization of the Effective Range 
Model using a momentum cutoff requires the constraint r s < 0. There is no 
apparent constraint on r s if we use dimensional regularization. This puzzle 



36 



was first pointed out in Ref. [IB] . The fact that dimensional regularization 
of the Effective Range Model gives the phase shift in Eq. (!85|) suggests that 
there might be a more elaborate model with a momentum cutoff that gives 
this phase shift. However, the results of Ref. [TB] imply that such a model 
cannot be a Hamiltonian quantum field theory. 



6.3 Time- derivative interactions 



Interactions that involve time derivatives are not allowed in a conventional 
Hamiltonian quantum field theory. However, they can be easily included in the 
Lagrangian approach. We will consider only the simplest Lagrangian model 
with a time-derivative interaction and show that it is equivalent to the Effec- 
tive Range Model. 

The Lagrangian density for a free atom is given in Eq. 0321) . The interaction 
term in the Lagrangian density for the Zero- Range Model is C in t = —Hint, 
where 7ii nt is given in Eq. (HoT) . The simplest generalization of the Zero-Range 
Model that includes a time- derivative interaction is a model whose interaction 
term is 



The combination of time and space derivatives acting on ip 2 is Galilean invari- 
ant. 

The diagrammatic analysis of the model whose Lagrangian is the sum of 
Eqs. ( 1321) and (1861) is identical to that in the Zero- Range Model except that the 
interaction vertex depends on the total energy E of the interacting atoms in 
their center-of-mass frame. The amplitude in this model can be obtained from 
the unrenromalized Green's function for the Zero-Range Model in Eq. ( 1511) 
simply by the substitution A — > A + AqK 2 , where k = {—mE — ie) 1 ^ 2 . The 
resulting expression for the Green's function is 



A(E) 



1 



An + XnK 2 47T 2 
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17) 



We obtain the T-matrix element T(k) by setting k = —ik. By using the 
expression for T(k) in Eq. (fT6j) and the expression for k cot 5o(k) in Eq. (EC 
we can determine the scattering length and effective range: 
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1 8-7T / 1 m A \ , 

- = — - + — A , 88a 
a m VAq 47T / 

167T An , , N 

r -—— f < 88b > 

By inverting these equations, we can determine how the bare parameters An 
and Ag must be tuned as functions of the cutoff A in order to keep a and r s 
fixed: 



A (A) = ( 1 - -aA] — , (89a) 

A / (A)=-r.^ri-^aA) 2 . (89b) 

Inserting these expressions into Eq. (J8"7|) and taking the limit A — > oo, we 
obtain the renormalized amplitude: 



(90) 

This is identical to the renormalized amplitude for the Effective Range Model 
in Eq. (ItjD"]) . Since all 2-body observables can be obtained from the Green's 
functions, we conclude that the Lagrangian field theory with the time-derivative 
interaction in Eq. fl86|) is equivalent to the Effective Range Model, at least in 
the 2-body sector. 
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6.4 Wilsonian renormalization 



The Wilsonian approach to renormalization requires finding a fixed-point 
Hamiltonian 7i* together with scaling operators (9* and their scaling dimen- 
sions 7„. The Hamiltonian can then be expanded around the fixed point as in 
Eq. (jl3"|) . The Wilsonian approach can also be applied to the larger space of 
Lagrangian quantum field theories. In this case, there is a similar expansion 
around a fixed-point Lagrangian £*. 

In a remarkable paper, Birse, McGovern, and Richardson carried out a com- 
plete Wilsonian analysis of 2-body S-wave interactions [20]. In a Lagrangian 
quantum field theory, the general S-wave 2-body operator can be expressed in 
terms of a potential V(p, k, k', A) that depends on the total energy E = p 2 /m, 
the squares k 2 and k 2 of the initial and final momenta, and the ultraviolet cut- 
off A. The requirement that the solution A(p, k, k') to the Lippmann-Schwinger 
integral equation be independent of the ultraviolet cutoff A gives renormaliza- 
tion group flow equations on the infinite-dimensional space of potentials. The 
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trivial fixed point is V = 0. There is also a nontrivial fixed point corresponding 
to the unitary limit of infinite scattering lengthQ 



V*(p,A) 
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p A + p 

A - o log 

2 A — p 



(91) 



This can be expanded as a power series in p 2 , so it corresponds to an interaction 
term with arbitrarily high-order time derivatives. A system infinitesimally 
close to the fixed point has a potential of the form V*(p, A) + v(p,k,k',A). 
The authors of Ref. [20] found all the scaling operators v(p,k,k',A). The 
scaling operators v (p, A) that depend only on the energy variable p have the 
scaling dimensions —1, 1, 3, 5, . . . that were deduced in Refs. [7f8l9] . There 
are also scaling operators v (p, k, k', A) that have nontrivial dependence on the 
momentum variables k and kl . They have scaling dimensions 2, 4, 5, 6, .... 
They may correspond to redundant operators in the Lagrangian formulation 
that have no counterparts in a Hamiltonian formulation and do not affect 
physical observables. 

Harada and Kubo and collaborators have applied the Wilsonian renormaliza- 
tion approach to Lagrangian field theories that are truncated in the derivative 
expansion [2~H22f23j . In Refs. [2TH22] they considered the 3-parameter La- 
grangian consisting of the free term in Eq. (l32l . the two interaction terms of 
the Effective Range Model in Eq. (1621) . and the time-derivative interaction in 
Eq. (1861) . This corresponds to a potential of the form 



V(p, k, k', A) = A (A) + Po (A)(A; 2 + k' 2 ) - X' (A)p 2 . (92) 

By demanding that the solution to the Lippmann-Schwinger integral equation 
have a cutoff-independent expansion to first order in p 2 , k 2 , and k 2 , they ob- 
tained renormalization group flow equations for the coupling constants Ao(A), 
Po(A), and X' . These equations have 3 fixed points: the trivial fixed point, the 
nontrivial fixed point, and an unphysical fixed point that gives complex scal- 
ing dimensions. Near the nontrivial fixed point, there is a relevant operator 
with scaling dimension —1 and two irrelevant operators with scaling dimen- 
sions 1 and 2. In Ref. [23J, the analysis was extended to the next order in 
the derivative expansion, where there are 8 coupling constants. They found 
three fixed points, two fixed lines, and a fixed surface. One of the fixed points 
was identified as the nontrivial fixed point associated with infinite scattering 
length. 



1 This differs from the expression in Ref. [2D] by a factor of 2 associated with 
identical particles. 
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7 Two-Channel Model 



A scattering model can describe multiple spin states of an atom. A minimal 
model has two spin states, including the one of primary interest. In general, 
there are 3 scattering channels corresponding to both atoms in the lower state, 
both atoms in the upper state, and atoms in both states. We consider the 
simple case of momentum-independent interactions between the atoms, as in 
the Zero-Range Model. We further simplify the problem by assuming that 
there is nontrivial scattering only in the two channels in which both atoms 
are in the same spin state. We will refer to this model as the Two- Channel 
Model. If the energies of atoms at rest in the two spin states are and z//2, the 
scattering thresholds are and v. We choose v > and we take the spin state 
of primary interest to be the one that is lower in energy. The parameters of 
the Two-Channel Model can be defined by specifying the S-wave phase shift 
in the lower spin channel to be 



k cot 5 (k) = 2~ 



1 



Ol2 



1/^22 + Vmu — k 2 



(93) 



In addition to the energy gap z/, there are three interaction parameters with 
dimensions of length: an, a±2, and 022- The interaction parameters are defined 
in such a way that the two channels decouple in the limit ai2 — > ±00. In this 
limit, an and a 2 2 reduce to the scattering lengths for the two independent 
channels. 



The scattering length and the effective range for the lower spin channel are 



J , (94a) 

2 

(94b) 

Note that the effective range r s is negative definite. There are various ways 
to tune the parameters to get a large scattering length a. If an < af 2 /a22, a 
large scattering length can be obtained by tuning the energy gap v to near 
the critical value (l/a22 — an/a^ 2 ) 2 /m. At v approaches the critical value, the 
effective range approaches r s — > o J 22o\ 2 l\ a \\{ a \\ a 22 — a ?2)]- A large scattering 
length can also be obtained by tuning the interaction parameter an to the 
critical value — a\ 2 \^/mv — 1/ 022] ■ 

A renormalizable Hamiltonian quantum field theory for a Two-Channel Model 
with two distinguishable particles with unequal masses was constructed by 
Cohen, Gelman, and van Kolck [24J . An essentially equivalent model has been 
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used to describe the effects of A A states on the two-nucleon system [25] . 
In Sec. I7.6[ we generalize the Two-Channel Model to one with iV scattering 
channels. 



7. 1 Hamiltonian 



The Two-Channel Model can be formulated as a quantum field theory with 
two complex fields i/j and ip2- The Hamiltonian density for the Two-Channel 
Model is the sum of a free term and an interaction term: 



Hint = t \ni^Hf + 7 V12 Wi>2? + (^W] + UoM4^) 2 - (95b) 



4 4 
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The Two-Channel Model is renormalizable in the 2-atom sector. The relations 
between the physical parameters and bare parameters Aq^- are 

1 Stt A , 22 2 a 

— = — t : rj— + -A, (96a) 

On ^ ^0,11^0,22 — An 12 7T 



v 0,12 

ai2 W Ao i llAo,22 — Aq j 12 

1 8tt A ,ii 2 



(96b) 

A. (96c) 



022 m Ao,nAo,22 — A i 2 vr 
These equations can be expressed more compactly in the form 



' ^(Aq 1 ) +-AS mn , (97) 

V / mn 7r 



dmn m K /mn 7T 

where Aq 1 is the inverse of the matrix Ao of coefficients Ao,^. 

Renormalization invariants are functions of the bare parameters that are equal 
to the same functions of the corresponding renormalized parameters indepen- 
dent of the value of the ultraviolet cutoff. Parameters that require no renor- 
malization, such as the energy gap u, are renormalization invariants. An al- 
ternative choice for the remaining parameters of the Two-Channel Model are 
the entries X mn of the matrix A whose inverse has the entries 

(A" 1 ) = j^-. (98) 

V / mn Kirn 
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With this choice of renormalized parameters, the traceless components of the 
matrices A -1 and Aq 1 are renormalization invariants: 

( A ") 12 =( A °V (99b) 

If we take the three independent interaction parameters to be these two renor- 
malization invariants and tr(Ag 1 ), the only bare parameter that must depend 
on the ultraviolet cutoff is t^Ag 1 ). 



7.2 Green's function 



The amputated connected Green's function A(E) for the Two-Channel Model 
is a 2 x 2 matrix that depends on the total energy E in the center-of-mass 
frame. The inverse of the matrix A(E) can be calculated analytically by solving 
the coupled system of integral equations represented diagramatically in Fig. HJ 



A(Eyl = ( -(V)ii " (m/47r 2 )(A - ttk/2) - (A ~ % \ 

{-(K 1 )^ - (Ao % - (m/47r 2 )(A - iry/mu + k 2 /2) ) ' 

(100) 

where k = {—mE — is) 1 / 2 . Inserting the relations in Eq. ( 1971) between the 
physical parameters and the bare parameters \o,ij, the dependence on the 
ultraviolet cutoff A disappears and the inverse matrix reduces to 



—1/an + K -l/ai2 
i —l/a\2 —l/a 2 2 + \rrrw + k 2 



(101) 



The square roots are defined for negative real arguments by the is prescription 
in the definition of k in Eq. (1521) . The explicit expression for the entries of this 
matrix are 



An{E) 
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h v mv + k z 
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+ K 

L an 


1 
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(102a) 
(102b) 
(102c) 
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These expressions agree with the results of Ref. 
associated with identical p articles q 



except for a factor of two 



7.3 T-matrix element 



The T-matrix element for the elastic scattering of atoms in the first channel 
with relative momentum k is obtained by evaluating An(E) in Eq. (I102ap at 
the total energy E = k 2 /m: 



_ /TN 8tt / 1 ., 1 
Mi(k) = — %k T 



h Vfriu — k 2 

a-22 



(103) 



Comparing with Eq. ffTrJj) , we obtain the equation for the phase shift in Eq. ( 1531 . 



7.4 Bound states 



The three entries of the matrix A(E) given by Eq. (11021) all have poles at the 
same values of the energy variable k. The values of k at the poles satisfy the 
equation 



(104) 

By solving for the square root and then squaring both sides of the equation, one 
can show that the roots of Eq. (11041) are also roots of a quartic polynomial. 
If the root Kj is real and positive, there is a stable bound state below the 
scattering threshold for the first channel with binding energy K 2 /m. If ^/mv > 
1/0-22, Eq. (11041) can have or 1 positive real roots. If a 2 2 > and y/mv < 
1/0-22, Eq. (11041) can have 1 or 2 positive real roots. Thus the number of stable 
diatomic molecules can be 0, 1, or 2. 

If the scattering length a is made large by tuning v to near its critical value, 
Eq. (11041) will have one small root k ~ 1/a. If a > 0, the corresponding bound 
state is the shallow dimer. The remaining roots of Eq. (11041) satisfy \k\ ^> l/N- 
As a — > ±00, these large roots approach solutions to a cubic equation: 



1 1 
— + — 



h \Jmv + 

O99 



k3 -^ 2 + + K + ^ [l- a f H (105) 

a ll \ a ll a 12 a 22 «12 / a 12 \a22 ai2/ 



2 There is an error in the component A22{E) in Ref. [24J, but it vanishes when the 
masses of the two particles are set equal. 
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Since this cubic equation was derived by squaring a square root, some of its 
roots may not correspond to roots of Eq. (11041) . If a n < 0, the limiting form of 
Eq. (11041) has no large positive real roots. If an > 0, it has one large positive 
real root. Thus the number of stable deeply-bound diatomic molecules is either 
or 1. 



7.5 Optical theorem 



If the energy E is real, the imaginary part of the amplitude in Eq. (I102al) is 




h Vitlu — k 2 — ie 

a-22 



9(E - v) 9(E) 



n -l 



S(E + K 2 Jm) 9{Ki 



(106) 

where the sum in the last term is over the 0, 1, or 2 positive real roots of 
Eq. (1 104ft . The expression for lm.An(E) is positive definite, in accord with the 
Optical Theorem in Eq. (|29|) . This is consistent with the absence of negative- 
norm states in the 2-atom sector of the Two-Channel Model. 



7.6 Multi-channel models 



The Two-Channel Model can be generalized to one with N spin states with 
nontrivial scattering only between pairs of atoms in the same spin state. We 
refer to this model as the N- Channel model. It has ^(N 2 + 3N — 2) parameters: 
the N — 1 energy gaps v n between the lowest 2-atom threshold iv\ = 0) and 
the n th 2-atom threshold and ^N(N + 1) interaction parameters a mn with 
dimensions of length. The parameters in the iV-Channel model can be defined 
by specifying the S-wave phase shift to be 



where C\ n (k 2 ) is the In cofactor of the N x N matrix whose mn entry is 
(m/87r)(— l/a mn + \/mv n — k 2 5 mn ). The scattering length and the effective 
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range are obtained from the low- momentum expansion of Eq. (11071) : 



\a u Gn(0) ~ ai n J 
n ^ C n (0)C[ n (0)- ^i(0)Cin(0) 

„=2 ainC/n^UJ 

where C[ n {k 2 ) = {d/dk 2 )C ln {k 2 ). 

The iV-Channel model can be formulated as a quantum field theory with N 
complex fields ip n , n = 1, • • ■ , N. The Hamiltonian density for the iV-Channel 
model is the sum of a free term and an interaction term: 



N r 1 1 



n=l L 



2m r " '' l 2 

N N 



(109a) 



-j 1\ 1\ 

Wint = tEE Ao, m „(^^n) 2 . (109b) 
m=l n=l 

Note that v\ = by definition. The iV-Channel model is renormalizable in 
the 2-atom sector. The energy gaps v n require no renormalization. The rela- 
tions between the physical interaction parameters a mn and the bare coupling 
constants A 0imn (A) are 

' HK (K 1 ) (HO) 



The amputated connected Green's function A(E) for this coupled-channel 
system is & N x N matrix that depends on the total energy E in the center- 
of-mass frame. The inverse matrix A(E)~ l can be obtained analytically by 
solving the coupled- channel integral equations represented diagramatically in 
Fig. HI After the renormalization using Eq. (IllOp . the mn entry of the inverse 
matrix reduces to 



(A(E)^ 



mn „ 

87T 



+ \ mv m + K 2 5 



where n = (—mE — ie) 1 ^ 2 . This inverse matrix can be inverted by using the 
identity 



^ - (n2) 
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where C nm {—n 2 ) is the nm cofactor of the matrix A(E)~ l given in Eq. fillip . 
Since the determinant is the sum of the products of the entries in any row (or 
column) and their cofactors, the mn entry of the matrix A(E) reduces to 



57T 



m 



C nm (mE) 



j^ C nl {mE) | 
1=1 



C nn (mE) 
C nm (mE) 



(113) 



The T-matrix element for elastic scattering in the first channel is obtained 
by evaluating the expression for An(E) given by Eq. (I113P at the energy 
E = k 2 /m: 



r xl {k) 



87T 

in 



ik — 



1 A C ln (k 2 ) 
Cn(k 2 )t 2 a ln 



-i 



(114) 



This T-matrix element gives the S-wave phase shift for the iV-Channel model 
in Eq. ffTOTD . 



8 Resonance Model 



The basic degrees of freedom in a scattering model can include molecular 
states as well as atoms. A minimal model has a single diatomic molecule 
that couples to a pair of atoms in the spin state of interest. We consider the 
simplest case of a momentum-independent interaction between the atoms and 
a momentum-independent coupling of the molecule to a pair of atoms. We 
treat the molecular degree of freedom as a point particle, so that the only 
structure of the molecule arises from its coupling to the atoms. We refer to 
this model as the Resonance Model. The parameters of the Resonance Model 
can be defined by specifying the S-wave phase shift of the atoms to be 



kcotSo(k) = --(\-^^-] . (115) 
m y mv — k z J 

The three parameters can be interpreted as the detuning energy v of the 
molecule relative to the 2-atom threshold, the coupling strength g of the 
molecule to a pair of atoms, and the self-coupling strength A of the atoms. 

The scattering length and the effective range are 
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m 2 is 2 



(116a) 



(116b) 



Note that the effective range is negative definite. If A = 0, the phase shift in 
Eq. (11151) has the same form as the phase shift for the Effective Range Model 
in Eq. ( !60l) . Thus the Effective Range Model with r s < is a special case of 
the Resonance Model. The scattering length in Eq. ( 1116a[) can be made large 
by tuning the energy gap v to near 0. The limiting value of the effective range 
is r s — ► — lQn/{m 2 g 2 ). 

The Resonance Model provides a natural description of atoms near a Feshbach 
resonance, where the dependence of the scattering length a on the magnetic 
field B has the form given in Eq. ( 1231) . That dependence can be reproduced by 
taking v to be linear in the magnetic field B while g 2 and A are proportional 

to ab g : 



v = -pL i {B-B i ), (117a) 

A = ^, (117b) 
m 

2 = _8 ! ra M 
m 

The parameter in Eq. (11171) can be interpreted as the difference between 
the magnetic moment of the molecule and twice the magnetic moment of an 
isolated atom. In general, both a^ g and can be slowly-varying functions of 
the magnetic field. 

The Resonance Model was first introduced by Kaplan as a model for non- 
relativistic particles with a large scattering length He used dimensional 
regularization which eliminated the need for explicit renormalization of the 
parameters. The Resonance Model was constructed independently by Kokkel- 
mans et al. [27] as a model for atoms near a Feshbach resonance. They derived 
the renormalization of the parameters that is required to make the observables 
independent of the ultraviolet momentum cutoff. In Sec. 18. 7\ we generalize the 
Resonance Model to one in which the atoms are coupled to iV different molec- 
ular states. 



8. 1 Hamiltonian 



The Resonance Model can be formulated as a quantum field theory with two 
complex fields: ip, which annihilates an atom, and 0, which annihilates a di- 
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atomic molecule. The Hamiltonian density for the Resonance Model is the 
sum of a free term and an interaction term: 



H free = — V# • W + ^V0 f ■ V0 + z/ o 0V, (H8a) 

H int = \go (0V 2 + ^ t2 0) + ^Ao(^V) 2 - (H8b) 

To avoid ultraviolet divergences, an ultraviolet cutoff A must be imposed on 
the momenta of virtual particles. 

The Resonance Model is renormalizable in the 2-atom sector, which consists 
of states containing two atoms or one diatomic molecule. The values for the 
bare parameters that are required to reproduce the phase shift in Eq. (11151) 
are 



A (A)=Z(A)A, 
g (A)=Z(A)g, 
u (A)=u-[l-Z(A)]g 2 /X 1 

where the renormalization constant Z is 



(119a) 
(119b) 
(119c) 



Z(A)=(l (120) 

Note that in the special case A = 0, which corresponds to the Effective Range 
Model, the only renormalization that is necessary is an additive renormaliza- 
tion of the detuning energy: 

u (A) = u + —g 2 A. (121) 

The bare parameter go is equal to its renormalized counterpart g. 

Using the renormalizations of the parameters given in Eqs. (I119p . one can 
construct two renormalization invariants: 



S/A = VA , (122a) 
u-g 2 /X = u -g 2 /X . (122b) 

If we take the three independent parameters to be the two renormalization in- 
variants in Eqs. (11221) along with A , the only parameter that must be adjusted 
as a function of the ultraviolet cutoff is Aq(A). 
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+ ■■■ 

Fig. 5. The series of diagrams that give the amputated connected Green's function 
for atom-atom scattering in the Resonance Model. 




Fig. 6. The integral equation for the amputated connected Green's function for 
atom-atom scattering in a model with a diatomic molecule that has a coupling to 
two atoms, such as the Resonance Model. 

8.2 Green's function 

The amputated connected Green's function for two atoms to evolve into two 
atoms can be calculated analytically in the Resonance Model by summing the 
series of diagrams shown in Fig. [5j Alternatively, it can be obtained by solving 
the integral equation illustrated in Fig. [61 This amplitude depends only on the 
total energy E of the two atoms in the center-of-mass frame: 



-i 

(123) 

where k = {—mE — ie) 1 ^ 2 . After inserting the renormalization conditions in 
Eqs. f)119p . the dependence on the ultraviolet cutoff A disappears and the 
amplitude reduces to 



A(E) 



A + 



E — Vq + is 
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8tt 



4tt 2 



A{E) 
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E — v + ie 
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8tt 



(124) 
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8.3 T-matrix element 



The T-matrix element T(k) for the scattering of atoms with relative momen- 
tum k can be obtained from the amplitude A(E) in Eq. (I124p by setting E 
equal to the total energy k 2 /m of the two atoms: 



T(k) 



2 \ -1 

mg z \ % 
-r H mk 

mv — k, 8n 



(125) 



Comparing with Eq. (1161) . we obtain the S-wave phase shift given in Eq. ( I115p . 



8.4 Bound states 



The amplitude A(E) in Eq. (I124p has poles in the energy variable k defined 
by Eq. (I521) . The values of n at the poles satisfy 



(x-^-X^-.. (126) 
\ mu + k 2 ) 8n K J 

This can be expressed as a cubic polynomial equation. The number of positive 
real roots is if v > and A < 0, 1 if v and A have the same sign, and 2 if 
v < and A > 0. Thus the number of stable diatomic molecules can be 0, 1 
or 2. 



If the scattering length a is made large by tuning v to near 0, Eq. (I126p will 
have one small root k ~ 1/a. If a > 0, the corresponding bound state is the 
shallow dimer. The two remaining roots of Eq. (I126p satisfy \k\ ^> l/|a|. As 
a — > ±00, these two large roots approach 



K± 
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mX 
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mXJ 
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X 



(127) 



If A < 0, neither of these roots is real and positive. If A > 0, the root k + is 
real and positive. It corresponds to a stable deeply-bound diatomic molecule. 



8.5 Optical Theorem 



If the energy E is real, the imaginary part of the amplitude in Eq. (1124=11 is 
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-1 

5(E + K 2 i /m)e{ Kl ), (128) 

where the sum in the last term is over the 0, 1, or 2 positive real roots Ki 
of Eq. 01261) . The expression for ImA(E) in Eq. f)128p is positive definite, 
in accord with the Optical Theorem in Eq. (T29]) . This is consistent with the 
absence of negative-norm states in the 2-atom sector of the Resonance Model. 
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8. 6 Resonance Model with ghost molecule 



The expression in Eq. (1116bj) for the effective range r s in the Resonance Model 
is negative definite. Thus this model cannot be used as a phenomenological 
description of atoms with positive effective range. Kaplan pointed out that 
this limitation can be avoided by taking the molecular field to be a ghost 
field whose quanta are states with negative norm [26]. The Lagrangian for this 
field theory differs from the Lagrangian for the Resonance Model only in the 
signs of a few terms. The Lagrangian density is the sum of a free term and an 
interaction term: 



•^free 
"^int 



' I- — 

. dt 2m 



a i/ 



\g (^ + ^) 



(129) 



where a = ±1. If a = +1, this is the Lagrangian density for the Resonance 
Model. If a — — 1, the free theory describes an atom and a ghost molecule. 
In the Feynman rules for the interacting theory, the only difference from the 
Resonance Model is in the propagator for the molecule. The propagator for a 
molecule with energy E and momentum P is i/[a(E — P 2 /(4m)) — v + ie]. 
The expression for the Green's function in Eq. (I124p is replaced by 



A(E) 



-i 



A + 
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a E — v + ie 
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-1 



(130) 



The S-wave phase shift in Eq. fll 15[) will be modified to 



*«**,(*) = -5: (a — ^yf 

m \ mv — a k z I 



(131) 
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The scattering length is given by Eq. (1116aj) for either sign a. The effective 
range differs from the expression in Eq. (1116bl) by a factor of a: 



T67T a X 




m?v 2 



(132) 



If a = —1, this is positive definite. 



The amplitude in Eq. (11301) has poles in the energy variable k defined by 
Eq. fl52|) . The values of k at the poles are roots of the equation 
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If a = — 1, the number of positive real roots can be or 2 if v and A have the 
same sign and 1 or 3 if they have opposite signs. If the energy E is real, the 
imaginary part of the amplitude in Eq. (I124j) is 



ImA(E) 




9(E) 

5{E + K 2 Jm)6{K i ), (134) 



where the sum in the last term is over the 0, 1, 2, or 3 positive real roots Hi 
of Eq. (11331) . Note that if a = —1, the last term need not be positive definite. 
Thus this model may have states that correspond to diatomic molecules with 
negative norm. 



The scattering length in Eq. (I116al) can be made large by tuning the energy 
gap v to be near 0. The limiting value of the effective range for a = — 1 is 
r s — > 167r/ (m 2 g 2 ), which is positive. There is one small root of Eq. (11331) that 
approaches k « 1/a in this limit. If a > 0, the corresponding bound state is 
the shallow dimer. For this root, the delta function contribution in Eq. (11341) 
has a positive coefficient, so the shallow dimer has positive norm. The two 
remaining roots of Eq. (11331) satisfy \k\ 3> As a — > ±oo, these two large 

roots approach 
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For the Resonance Model with a ghost molecule, we set a = —1. A root n 
that is real and positive corresponds to a stable deeply-bound state. If the 
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root is in the range < k < m 2 g 2 / (Art) , the delta function contribution in 
Eq. (11341) has a negative coefficient, so the bound state is a negative-norm 
state. If A < 0, the only large positive real root is k + , and it corresponds to 
a negative-norm deeply-bound bound state. If < A < 16n 2 / (m 3 g 2 ) , both 
n + and K- are positive real roots. They correspond to a positive-norm and a 
negative-norm deeply-bound state, respectively. If A > IQn 2 / (m 3 g 2 ) , there are 
no large positive real roots, so there are no stable deeply-bound states. Even 
if it includes deeply-bound negative-norm states, the Resonance Model with 
a = — 1 may still be useful as an approximate description of atoms with a large 
scattering length and a positive effective range at energies small compared to 
1/ (mr 2 ). 



8.7 Multi-resonance models 



The Resonance Model can be generalized to one with N diatomic molecules 
that couple to a pair of atoms but not to each other. We refer to this model as 
the N -Resonance Model. It has 2iV-|-l parameters: the self-coupling strength A 
for the atoms, a detuning energy v n for each of the iV molecules, and a coupling 
strength g n for each of the N molecules. The parameters in the iV-Resonance 
Model can be defined by specifying the S-wave phase shift to be 

kcat5 (k) = -—(\-j2 m9l 7 2 ) . (136) 

m \ n^l mz/ « - « / 

The scattering length and the effective range are 
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Note that the effective range is negative definite. 

The Hamiltonian density for the iV-Resonance Model is the sum of a free term 
and an interaction term: 
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As we will show below, the relation between the bare coupling constants 
in Eq. (11381) and the renormalized parameters in Eq. (11361) can be expressed 
compactly as an equality between two functions of the energy E: 



(139) 



By taking the limit E — > oo, we find that the relation between Ao and A has the 
same form in Eq. (I119al) as in the Resonance Model. Both sides of Eq. (11391) 
can be expressed as ratios of polynomials in E of degree N. The equality 
between these rational functions can be expressed as 2N coupled nonlinear 
equations for the remaining 2N parameters of the A^-Resonance Model. 

The amputated connected Green's function for two atoms to evolve into two 
atoms can be calculated analytically by solving the integral equation repre- 
sented diagramatically in Fig. The amplitude A depends only on the total 
energy E in the center-of-mass frame: 
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where k = {—raE — ie) 1 / 2 . Inside the square brackets in Eq. (11401) . the sum 
of the term raised to the power —1 and the term (m/47r 2 )A can be expressed 
as the ratio of two ATth order polynomials in E. By the fundamental theorem 
of arithmetic, the polynomial in the numerator has A" zeroes in the complex 



plane, which we denote by v n 



n 



N . Except on sets of the bare param- 



eters that have measure zero, those zeroes will be distinct. By decomposing 
the inverse of that ratio of polynomials into partial fractions, the amplitude 
in Eq. (11401) can be expressed in the form 



A(E) 



A + 



N 

E 



-i 



-i 



Z^E - v n + ie i 



m 
8tt' 



-K 



(141) 



By comparing Eqs. (11401) and f 1 1 4 1 j) . we obtain the renormalization condition 
in Eq. fTl39l . 

The condition that the Hamiltonian be hermitian requires that the bare pa- 
rameters Ao, vo,n and go,n i n the Hamiltonian density in Eqs. (11381) be real 
valued. The expression in Eq. (I140p is therefore a real- valued function for real 
K. This implies that the parameters v n and g 2 in Eq. (11411) are either real or 
they come in complex conjugate pairs. This is sufficient to guarantee that the 
phase shift in Eq. (I136p is real valued, which is necessary if the model is to 
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describe atoms that have no inelastic scattering channels. To show that the 
iV-Resonance Model is renormalizable in the 2-atom sector for arbitrary real 
values of the parameters v n and g n , one must show that the renormalization 
conditions in Eq. (I139p have real solutions for the bare parameters i/ , n an d 
<7o,n as A ^ oo. We will verify this below for the Two-Resonance Model. We 
have not shown that this is the case for the iV-Resonance Model with N > 3. 

The renormalization conditions for the Two-Resonance Model are obtained 
by setting N = 2 in Eq. (11391) . We can use Eq. (j!20p to express the ultraviolet 
cutoff A SIS db function of A and Z: 




(142) 



If we also use Eq. (I119al) to eliminate Ao in favor of A and Z, the renormaliza- 
tion condition in Eq. (11391) reduces to 



(£-z/i)(£-i/ 2 ) 

(E - vi){E - v 2 ) + [gl(E - v 2 ) + g 2 2 (E - u x )]/X 

_ (E - u 0A )(E - upj) + (Z- 1 - Z-^jgl^E - u^) + gl 2 {E - u 0A )]/\ 
(E- u Q ,i)(E - u 0>2 ) + Z-^gl^E - z/ , 2 ) + gl 2 (E - i/ 0j1 )]/A 

(143) 

We obtain four coupled equations by matching the constant terms and the 
linear terms in E in both the numerator and denominator: 



v\v% = ^0,1^0,2 - Z Z l ){gl jX v Gj2 + ^0,2^0,0/A, (144a) 

vx + V2 = v ,i + "o,2 - Z-\l - Z- l \g\ x + gl 2 )/\, (144b) 

"1^2 ~ (g\v2 + g\v\)l\ = ^0,1^0,2 - Z~ l (^0,1^0,2 + ^O^O.O/A, (144c) 

v\ + ^2 - {g\ + gl)/X = vo,i + "o,2 ~ Z-\gl x + gl 2 )/\. (144d) 

These equations can be solved for the bare parameters as functions of the 
renormalized parameters and Z: 



^o,i = 2 ("i + ^2 + C) 

"0,2 = ~ (^1 + ^2 " C) 



9o,i 

90,2 



Z ( „2 , 2 



^ ( Jl , „2 



1 


(145a) 


? 


(145b) 


mi/a 
c) ' 


(145c) 


D\ 1/2 

c) ' 


(145d) 
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where 



v x = v 1 -(\-Z)g\j\ (146a) 
v 2 = V2 - (1 - z) g\j\ (146b) 

C = [{v x - v 2 f + 4(1 - Zf glgi/X 2 ] 1/2 , (146c) 

D={v 1 - V2M - gl) - 4(1 - Z) glgl/\. (146d) 

Note that the expressions for the bare parameters f ,i, g ,i> u 0>2 , and g ,2 m 
Eq. (I145p are each separately symmetric under the interchange of the renor- 
malized parameters: V\ <-> v 2 and g\ «-> (72- The Two- Resonance Model is 
renormalizable for given values of 1/1, gi, u 2 , and #2 if the solutions for the 
bare parameters in Eq. f)145|) are real valued in the limit A — > 00. Since the 
expression for C in Eq. (1146cj) is the square root of a manifestly positive quan- 
tity, it is clear that the parameters 1/0,1 and ^0,2 in Eqs. (I145aj) and fl!45b|) are 



real. The condition that the parameters #0,1 and c/0,2 in Eqs. (j 145 eft and (I145dl) 
are real is g\ + g\ > \D\/C. It is easy to verify that this is satisfied for any 
real values of u±, u 2 , g%, and g 2 , even at finite A. 



The Two- Resonance Model has been previously discussed in Ref. [27J, but 
their results for the renormalization of the parameters of the model are incor- 
rect. They obtained the correct result in Eq. (1119 al) for the renormalization of 
A. However, their results for the renormalizations of the parameters v\ and g\ 
for the first resonance are given by the same simple equations in Eqs. (1119bl) 
and fll 19c|) as in the Resonance Model, while their results for the renormaliza- 
tions of the parameters v 2 and g 2 for the second resonance are given by more 
complicated equations. It is easy to see from symmetry considerations that 
these results cannot be correct. The Hamiltonian of the N- Resonance Model 
is symmetric under the interchange <pi <-> <p 2 of the two molecular fields if we 
also interchange the bare parameters: z/0,1 ^ ^0,2 and go,i ^ #0,2- Thus the 
expressions for the bare parameters v\q and gi t o should differ from those for 
z/ 2) o and g 2 $ simply by interchange of the renormalized parameters. Our ex- 
pressions for the renormalized parameters in Eqs. (I145p satisfy this condition, 
but those in Ref. [271 do not. 



The amplitude in Eq. (11411) has poles in the energy variable k defined by 
Eq. (|52|) . The values of k at the poles are roots of the equation 



TV 2 \ _1 

mg n \ m 



_^z^ = ( 14 7) 
If the energy E is real, the imaginary part of the amplitude in Eq. (I14ip is 
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ImA(E) 



mk 



( 



n=l 



N 



mvrn - k 2 





4ir(mv n + k 2 ) 2 




5(E + K 2 Jm) 9{Ki) 



(148) 



where the sum in the last term is over the positive real roots K{ of Eq. (j!47p . 
This expression is positive definite, in accord with the Optical Theorem in 
Eq. fl29|) . This is consistent with the absence of negative- norm states in the 
2-atom sector of the iV-Resonance Model. 



Acknowledgments 

This research was supported in part by DOE Grants DE-FG02-91-ER4069, 
DE-FG02-05ER15715, and DE-FG02-06ER41449. 



A Integrals 

Some of the results in this paper require the evaluation of integrals over the 
momentum k of virtual atoms of the form 



where n is an integer. These integrals are functions of the total energy E of a 
pair of atoms in the center-of-mass frame. It is convenient to express them in 
terms of the variable k = (—mE — ie) 1 / 2 defined in Eq. (1521) . For n > 0, the 
integrals given by Eq. flA.lj) are ultraviolet divergent. They can be regularized 
by imposing a momentum cutoff: < A. The integrals then satisfy a simple 
recursion relation: 



This can be used to express the integrals l2 n {E) for n > 1 in terms of Io(E), 
whose value is 




(A.l) 



2tt 2 2n + 1 



n 2 I 2n _ 2 {E). 



(A.2) 




(A.3) 
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The expressions for the integrals are particularly simple if we take A to be so 
much larger than [m\E\) 1 ^ 2 that we can neglect terms that decrease as inverse 
powers of A/ '{m\E\) l l 2 as A — > oo. The integral in Eq. (1A.3j) then reduces to 

WB) = _| ; (A-| K ). (A.4) 
The next two integrals in the sequence are 



(A.5a) 
(A.5b) 



With dimensional regularization, the integral in Eq. (1A.1I) is generalized to an 
integral over a space with dimension 3 — 2e and then analytically continued 
to e = 0. This procedure automatically subtracts any power ultraviolet diver- 
gences. The resulting expressions for the integrals can be obtained simply by 
setting A = in Eqs. OA. 4ft and (IA.5D : 



Tfl 

I (E) = —k, (A.6a) 

I 2 (E) = --k\ (A.6b) 

Tfi 

I 4 (E) = — k . (A.6c) 

4:71 
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